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Chsmical  Profiling  of  Silicon  Nitride  Structures 


Abstract 

X-ray  phoioeiectron  spectroscopy  (XPS),  spectroscopic  ellipsometry 
(SEX  and  scanning  electron  microscopy  (SEM)  have  been  used  to  study 
structural  and  chemical  inhomogeneities  in  several  electronic  materials  and 
device  structures  of  relevance  to  radiation  hard  electronics.  The  systems 
studied  include  metal-nitride-oxide-semiconductor  (MNOS)  structures,  silicon 
oxynitride  (SiOxNy)  formed  by  the  thermal  nitridation  of  S1O2.  and 
semiconductor-on-insulator  (SOI)  structures.  Studies  of  MNOS  structures 
suggest  that  the  effect  of  H2  annealing  is  to  make  the  Si3N4/Si02  interface 
less  abrupt  by  causing  interdiffusion  of  silanol  and  silamine  groups  with 
subsequent  oxynitride  formation.  Another  effect  of  the  annealing  appears  to 
be  to  relieve  the  strain  at  the  Si02/Si  interface. 

The  chemical  and  structural  inhomogeneities  in  SiOxNy  formed  by  the 
thermal  nitridation  in  NH3  of  SiC^/Si  have  been  studied.  The  nitrogen 
concentration  depth  profile  as  a  function  of  nitridation  time  and  temperature 
has  been  measured  with  XPS.  The  peak  interfacial  nitrogen  concentration 
initially  occurs  at  the  SiOxNy/Si  interface,  then  appears  to  move  away  from 
the  interface  with  increasing  nitridation  time  for  sufficiently  high  nitridation 
temperatures.  The  kinetics  of  nitridation  has  been  modeled,  and  reasonable 
agreement  is  obtained  with  the  experimental  results,  provided  the  effects  of 
interface!  strain  are  explicitly  taken  into  account.  Electrical  data  in  the 
literature  are  found  to  be  correlated  with  the  nitrogen  distributions  measured 


iii 


here,  and  are  explained  within  a  model  based  on  the  strain-dependent 
formation  of  defect  states. 

Studies  of  SO!  structures  formed  by  ion  implantation  show  that  the 
cuned  oxide/Si  interfaces  are  smoother  with  mgner  implantation 
temperatures  in  the  range  400'  C  to  525'  C.  However,  the  stoichiometry  of 
the  native  oxides  and  thermal  SiC^/Si  interface  suggests  that  a  higher 
implantation  temperature  results  in  a  rougher  surface,  which  may  possibly 
affect  device  characteristics.  SOI  structures  formed  by  zone-melting 
recrystailization  have  also  been  studied,  but  these  studies  are  incomplete  at 
this  time. 
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2. 1 .  Schematic  diagram  of  the  system  used  for  gas  phase  XeF2 

etching  of  Si . 

3.1  NIc  and  0  1  s  photoeiectrcn  peak  intensities.  IW1  s  and 

l0ls,  measured  during  the  course  of  depth  profiling  sample 
corrected  for  photoionization  cross  section  differences . 

3.2.  IN1S/(I01S  +  IN1  '■  for  (a)sample  1.  and  (b)  sample  2, 

correaed  for  pnotoionization  cross  section  differences . 

3.3.  Distribution  of  N  in  SiO  N  produced  from  a  ICO  A-thick 

SiO,  film  grown  in  dry  02at  1000'C,  then  nitrided  at 
1000‘C  (a)  10  minutes,  (b)  30  minutes,  (c)  120  minutes.  . 

3.4.  Distribution  of  N  in  SiOJM  produced  by  nitriding  a  ICO 
A-thick  Si02  film  at  1 1 50‘C  for  (a)  5  minutes,  and 

(b)  10  minutes . 

3.5.  (a)  Comparison  of  the  distribution  of  N  in  SiOxN  films 
produced  from  100  A-thick  SiO»  films  grown  in  ary  02  at 
825'C  (solid  circles)  and  1000'C  (open  circles),  then 
nitrided  at  800‘C  for  60  minutes,  (b)  Distribution  of  N 

in  a  SiOxN  film  produced  from  a  100  A-thick  Si02 

film  nitrided  at  800‘C  for  240  minutes . 

3.6.  (a)  Magnitude  of  the  fiatband  voltage  shift  |AVFBj  as  a 

function  of  mtridation  time  for  600-A-thick  oxides  nitrided 
at  1000'C  (from  Ref.  57).  (d)  Peak  interfacial  nitrogen 
concentration  [N]  and  distance  of  [N]  from  the  interface 
as  functions  of  tN,  for  a  100-A-thick  oxide  nitrided  at 
1000'C  (from  Fig.  3.3) . 

3.7.  Schematic  representation  of  the  distribution  of  defects 
at  various  nitridation  times. 

Inset:  Schematic  diagram  of  the  free  energy  of 
formation  AGF  vs.  configuration  coordinate  for  the 
formation  of  SiO  N  /Si  without  defects  (solid  line)  and 
with  defects  (dashed  line)  to  relieve  strain  energy . 

3.3.  Surface  fluorine  concentrations  measured  during  the 
course  of  depth  profiles  of  samples  nitrided  at  1000'C 
for  (a)  10  minutes,  and  (b)  30  minutes . 
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3.9.  Surface  fluorine  concentrations  measured  during  the 
course  of  depth  profiles  of  samples  nitrided  at  1 15 O'C 

for  (a)  5  minutes,  and  (b)  10  minutes . 

3.10.  Surface  fluorine  concentrations  measured  during  the 
course  of  depth  profiles  of  samples  (a)  oxidized  at  825‘C, 
ther  mtnoed  a;  800‘C  for  60  minutes,  (b)  oxidized  at 
iGOC'C,  tnen  nitrided  at  800'C  for  60  minutes,  and  (c) 
oxidized  at  1000‘C,  then  nitrided  at  SOO'C  for  240  minutes.  ... 

3.11.  Schematic  overview  of  the  kinetic  model . 

3  -  2 .  Calculated  distribution  of  N  in  a  SiOxNv  film  produced 

by  nitriding  a  92  A-thick  Si02  film  at  1000‘C  for  (a)  5  minutes, 

(b)  10  minutes,  (c)  20  minutes,  and  (d)  30  minutes . 

3.13.  (a)  Calculated  curve  from  Fig.  3.12(d). 

(b)  Effect  of  including  an  activation  energy  for  nitrogen 
incorporation  in  the  oxide  which  decreases  linearly  within 
26  A  of  the  interface  as  the  interface  is  approached. 

(c)  Same  as  (b),  but  including  activation  energies  for 
hopping  of  diffusing  species  in  the  oxynitride  which 
decrease  linearly  in  the  interfacial  region. 

(d)  Same  as  (c)  but  extending  the  region  of  graded 
reactivity  to  include  the  two  Si  layers  in  the  substrate 

which  are  closest  to  the  SiOxNy/Si  interface . 

3.14.  Comparison  of  the  experimental  nitrogen  distribution 
m  a  SiOxN  film,  produced  from  a  100  A-thick  SiO,  film 
nitrided  tor  30  minutes  at  1000’C  (squares,  from  Fig. 

3.3(b)),  to  the  calculated  nitrogen  distributions  (a)  from 

i  iy.  3.13(et/,  anu  from  rig.  C.  i  3(d) . 

3.15.  Comoarison  of  the  experimental  nitrogen  distribution 
in  a  SiO  N  film,  produced  from  a  100  A-thick  Si02 
film  nitrided  for  10  minutes  at  1150‘C  (squares,  from 
F'g  3  4(b)1,  to  the  calculated  nitrogen  distribution 

(a)  using  the  same  kinetic  parameters  ac  in  Fig.  3.12(d), 
and  (b)  including  the  effects  of  a  reversible  mtridation . 

3  16.  Comparison  of  the  experimental  data  points  from  Fig. 

3.3(b)  (squares)  to  the  calculated  nitrogen  distribution 

using  the  same  kinetic  parameters  as  in  Fig.  3.15(b) 

for  a  30  minute  mtridation  at  1000'C . 


3.  ‘  ~  Ccrroanson  of  the  experimental  nitrogen  distribution  m 
SiOxNv  fiims,  produced  from  100  A-thicK  SiO-  films 
oitriaecj  for  60  minutes  at  800‘C  (solid  and  open  circles, 
from  Fig.  3.5(a)),  to  the  calculated  nitrogen  oistnbution 
(.a)  using  kinetic  parameters  appropriate  for  low 
temperature  mtriaation,  anc  (b)  using  kinetic  parameters 
aopropriate  for  high  temperature  mtndation . 

3  ’8.  Calculated  nitrogen  distribution  using  the  same  kinetic 
carameters  as  Fig.  3.13(dl  for  a  30  minute  mtndaticn 
at  1000'C  for  an  initial  oxide  thickness  of  ( a )  100  A, 

(b)  1 50  A,  and  (c)  200  A . 

3  "3  Soect.a  of  tan  mj  obtained  for  (a)  a  1030  A-thick  layer  of 
SiO-,  on  Si.  (bi  a  1030  A-thickiayer  o*  SO?  on  Si. 
nitriaed  for  2  nours  at  800‘C,  and  icl  a  1030  A-thick 
layer  of  Si02  on  Si.  mtnoed  for  2  hours  at  1 000'C . 

3.20.  Spectra  of  cos  A  corresponding  to  the  data  in  Fig.  3  19 . 

3.2'  Best  fit  to  the  data  from  Figs.  3.19(c)  and  3.20(c),  obtained 
using  a  single-layer  model  in  which  the  dielectric  consists 
of  a  physical  mixture  of  Si02  and  Si3N4.  Only  the  fit  to 
cos  A  is  shewn,  since  tan  u/  could  be  fit  well  enough  that  no 
difference  between  the  data  and  the  fit  is  apparent  to  the  eye. 

3.22  Comparison  of  tan  \j  spectra  obtained  from  sample  1 1 

using  aperture  sizes  of  0.6  mm  and  2  mm . 


3.23  Comparison  of  cos  A  spectra  obtained  from  sample  1 1 

using  aperture  sizes  of  0.6  mm  and  2  mm .  82 

3.24.  Calculated  spectrum  of  tan  v  for  a  1030  A-thick  dielectric 
on  Si,  consisting  of  a  physical  mixture  of  10%  Si3N4  and 

90%  SiQ2,  in  the  energy  range  0.5-5. 7  eV .  34 

3.25.  Calculated  spectrum  of  cos  A  for  a  1030  A-thick  dielectric 
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4.1 .  SEM  micrographs  of  the  surface^  of  SOI  samples  formed 
by  ion  implantation  of  oxygen  at  (a)  5C0CC,  ncanneal. 

(b)  425=C,  2  hour  1150°C  post-implantation  anneal, 

i c)  500°C,  2  hour  1 150°C  anneal,  and  (d)  525=0,  2  hour 
1 150=C  anneal . 

4.2.  XPS  SDectra  of  the  Si  2 p  region  for  the  samples 

corresponding  to  the  data  in  Fig.  4. 1 . 
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SEM  micrograons  (magnification  10-1 3kx)  of  the  surfaces 
o?  SOI  samples  formea  by  ion  .replantation  of  oxygen  at 
(a:  500;C.  no  anneai,  (b)  425°C,  2  hour  1150:C 
ocst-imolantation  anneal,  (c)  500°C,  2  hour  i150°C  anneai, 
and  id)  525'C.  2  hour  1150°C  anneal.  The  surface  Si  iayer 
nas  oeen  removed  by  XeF2  etching,  exposing  the  buried  oxide 

Same  as  Fig.  4.3,  but  with  a  magnification  of  65kx . 

SEM  micrographs  (magnification  50kx)  of  tne  surfaces  of 
SOI  samples  formed  by  ion  implantation  of  oxygen  at 
(a)  500°C,  no  anneal,  (b)  4252C,  2  hour  1 1 50°C 
post-impiantation  anneal,  (c)  500°C,  2  hour  1150CC  anneal, 
and  id)  525°C,  2  hour  1150°C  anneal.  Tnfe  surface  Si  layer 
nas  been  removed  by  XeF?  etching,  and  the  buned  oxicc 
layer  has  been  removed  wlTh  HF,  exposing  the  buned 
oxide'substrate  interface . 

Same  as  Fig.  4.5,  but  with  a  magnification  of  200kx . 
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Chapter  1 

introduction 

Tne  resu:;s  of  studies  of  matenais  and  device  structures  of  reievance  tc 
raciation  nard  e'ectronics  is  reported  herein.  This  work  is  a  follow-on  to  a 
research  project  VJPL  Task  Plan  #90-1276)  performed  by  JPL  for  RADC  anc 
involves  a  coordinated  study  aimed  at  understanding  the  chemical  anc 
structural  features  which  limit  the  radiation  stability  of  devices,  oamculariy 
those  based  on  silicon  nitride  and  oxynitride  dielectric  films.  The  participants 
have  ioeen  Westinghouse,  JnL,  MIT  Lincoln  Laboratory,  and  the  University  of 
Surrey,  England,  UK.  The  structures  studied  include  metal-nitride-oxioe- 
semiconductor  (MNOS)  structures,  silicon  oxynitride  (SiOxNy)  films  formed  by 
the  thermal  nitridation  of  SiOg  films  in  NH3,  and  semiconductor-on-insulator 
(SOI)  structures.  The  objectives  of  the  study  are  to  determine  the  chemical 
structure  of  silicon  nitride  and  oxynitnde  and  their  interfaces  with  Si  and  S iOg, 
to  identify  the  key  defects,  and  to  model  the  mechanism  of  defect  generation 
which  may  limit  device  penormance. 

Many  proposed  electronic  device  structures  require  insulator  or 
semiconductor  thin  films  which  are  only  a  few  lattice  constants  thick.  For 
example,  megabit  dynamic  random  access  memory  devices  require  an 
insulator  thickness  near  100  A.  In  the  past,  defects  and  inhomogeneities  on 
an  atomic  scale  could  often  be  ignored  since  they  were  so  much  smaller  than 
device  dimensions.  As  device  dimensions  become  smaller,  an  under¬ 
standing  of  interfacial  prooerties,  defects,  and  inhomogeneities  on  an  atomic 


scale  becorres  increasingly  important. 

Defects  and  inhomogeneities  can  be  classified  as  being  either  structural 
or  cnemica..  Structural  defects  include  extenaed  defects  such  as  gram 
ocuncaries  and  dislocations,  anc  defects  on  an  atomic  scale,  sucn  as  bond 
angle  distortions,  voids,  and  dangling  orbitais.  Chemical  inhomoaeneities 
include  inhomogeneous  distributions  of  constituents  in  a  random  alloy  (alloy 
clustering),  a  logical  extension  of  which  is  phase  separation  in  the  context  of 
metastabie  alleys.  In  this  work,  neitner  extended  structural  defects  nor 
chemical  inhomogeneities  such  as  trace  impurities  will  be  considered. 

It  is  to  be  noted  that  any  interface  between  two  different  materials  is 
mnerently  inhomogeneous,  both  structurally  and  chemically.  Structural 
inhomogeneities  arise  from  lattice  mismatch  between  the  two  materials, 
which  can  be  accomodated  by  lattice  distortion  of  one  or  both  materials,  or  by 
defects  which  allow  Doth  materials  to  retain  their  preferred  lattice  dimensions. 
Chemical  inhomogeneities  arise  from  the  fact  that  the  interface  is  a  boundary 
between  chemically  distinct  phases,  so  that  interfacial  atoms  necessarily 
have  bonding  which  is  different  than  that  found  in  either  bulk  material. 

Amorphous  materials  are  inherently  structurally  inhomogeneous,  in  the 
sense  that  the  atoms  do  not  occupy  well-defined  lattice  sites,  i.e.  the  material 
is  topologically  disordered.  Amorphous  materials  can  be  chemically 
homogeneous,  however.  A  chemically  homogeneous  amorphous  material 
which  has  been  extensively  studied''3  is  SiC^,  which  consists  of  Si 
fetrahedra  connected  by  bridging  oxygens  in  a  ring  structure.  Bulk 
amorohous  S1O2  (a-SiC^)  consists  primarily  of  6-membered  rings, 


corresoondmg  to  an  Si-O-Si  bridging  bond  angle  of  144\  the  same  as  that 
found  in  a-quariz.  Amorphous  SiOg  differs  from  a-quartz  in  having  a 
significant  void  fraction,  manifested  in  a  lower  density  (2.65  g/cm3  for  a-quartz 
vs.  2.2  g.c.m2  for  a-SiOg},  in  a  distribution  of  Si-O-Si  bond  angles,  and  hence 
in  a  distribution  of  ring  sizes.  These  distributions  are  affected  by  strain: 
soecificaily,  the  strain  which  is  localized  within  30  A  of  the  SiOg/Si  interface  in 
thermally  grown  oxides  on  Si  has  been  shown'-2  to  result  in  a  higher 
proportion  of  4-memoered  rings  fSi-O-S.  bridging  bone  angle  of  120”;  m  me 
interfaciai  region  compared  to  the  unstrained  bulk  SiOg.  The  results  of  these 
studies1-3  have  important  imolications  for  this  work,  as  will  be  seen. 

Another  amorphous  material  which  is  ideally  homogeneous  chemically 
is  Si3N4.  The  study  of  CVD  Si3N4  and  its  interface  with  SiOg  in  MNOS 
structures  was  the  primary  work  of  the  previous  task,  and  the  results  are 
described  elsewhere.4'9  Some  follow-on  work  has  been  done  in  this  task, 
concerning  the  effects  of  Hg  annealing  on  the  S i 3 N 4  and  the  SigN^SiO? 
interface,  and  the  results  are  described  herein. 

Amorphous  alloys  may  exhibit  chemical  inhomogeneities  in  addition  to 
the  structural  inhomogeneities  which  are  inherent  in  amorphous  materials.  In 
this  work,  SiOxNy  is  an  example  of  a  material  with  both  chemical  and 
structural  inhomogeneities.  The  SiOxNy  has  been  formed  by  the  thermal 
nitridation  in  NH3  of  SiOg  films  on  Si.  As  previously  mentioned,  the 
structures  of  SiOg  and  the  SiOg/Si  interface  have  been  extensively  studied, 
so  that  the  initial  state  of  the  sample  prior  to  nitridation  -has  been 
characterized.  Any  chemical  inhomogeneities  which  arise  during  the 


mtndation  can  be  followed  by  measuring  the  nitrogen  distribution  in  the  film  at 
various  nitndation  times.  It  may  thus  be  possible  to  elucidate  the  mechanism 
and  kinetics  of  tne  nitndation  process,  and  to  relate  the  chemical 
mnomogeneities  to  tne  Known  structural  inhomogeneities  present  in  the  initial 
S1O2  film.  Specifically,  the  strain  which  is  present  at  the  SiC^/Si  interface 
may  affect  the  nitridation  kinetics.  These  studies  have  been  done,10'18  and 
the  results  are  described  here. 

Two  types  of  SOI  structures  have  aiso  been  studied.  The  first  type  is 
formed  by  the  high  dosage  implantation  of  O  ions  into  a  crystalline  Si  wafer  at 
eievated  temperatures,  followed  by  a  high  temperature  anneal  to  redistribute 
the  implanted  oxygen  and  form  a  buried  SiC>2  layer.  Samples  of  this  type 
were  obtained  from  Dr.  Peter  Hemment  of  the  University  of  Surrey,  England. 
Studies  on  similar  samples  have  been  done,  and  the  results  are  reported 
elsewhere.19'22  The  second  type  of  SOI  structure  studied  is  formed  by 
growing  a  thick  thermal  Si02  layer  on  a  crystalline  Si  substrate,  depositing 
polycrystalline  Si  on  top,  and  depositing  a  thick  Si02  capping  layer  over  the 
poly-Si.  The  entire  structure  is  heated,  and  a  graphite  strip  heater  is  passed 
over  the  top  to  melt  the  poly-Si  and  recrystallize  it  as  a  single  crystal  Si  layer. 
Samples  of  this  type  were  obtained  from  Dr.  C.  K.  Chen  of  MIT  Lincoln 
Laboratory.  Studies  on  similar  samples  have  been  done,  and  the  results  are 
reported  elsewhere.23 

In  order  to  study  the  problems  discussed  above,  three  techniques,  which 
are  in  many  ways  complementary,  have  been  utilized:  x-ray  photoelectron 
soectroscoDy  (XPS),  scanning  electron  microscopy  (SEM),  and  spectroscopic 
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eiliDSometrv  (SE).  SEM  is  a  standard  technique  for  imaging  morphological 
features  on  surfaces.  A  detailed  description  of  the  technique  and  its 
applications  can  be  found  in  many  texts.24  and  will  not  be  described  here. 
XPS,  also  Known  as  ESCA  (electron  spectroscopy  for  chemical  analysis),  is 
based  on  the  photoelectric  effect,  in  which  photons  incident  on  a  specimen 
are  absorbed  and  electrons  are  emitted.  Only  a  brief  description  of  the  basic 
concepts  is  presented  below,  since  excellent  reviews  are  available  in  the 
literature.25'23  The  XPS  spectrometer  used  in  this  work  is  describee 
elsewhere. 1 

The  binding  energy  of  an  electron  relative  to  the  vacuum  level  is  the 
energy  required  to  remove  the  electron  to  infinity  with  zero  kinetic  energy. 
The  measured  kinetic  energy  of  the  emitted  electron  is  thus  equal  to  the  x-ray 
photon  energy  minus  the  binding  energy  of  the  state  from  which  the  electron 
was  emitted.  Measurement  of  the  kinetic  energy  of  the  emitted  electrons  will 
thus,  in  principle,  directly  yield  the  binding  energies  of  the  occupied  states.  In 
practice,  binding  energies  in  XPS  are  measured  relative  to  the  Fermi  level, 
rather  than  the  vacuum  level.  In  this  case,  the  kinetic  energy  of  the  emitted 
electron  is  equal  to  the  x-ray  photon  energy  minus  the  sum  of  the  binding 
energy  relative  to  the  Fermi  level  and  the  work  function  of  the  spectrometer. 
Each  element  has  a  set  of  energy  levels  (corresponding  to  atomic  and 
molecular  orbitals)  with  binding  energies  characteristic  of  that  element.  In  a 
photoelectron  intensity  vs.  binding  energy  plot,  a  set  of  peaks  will  be 
observed  which  correspond  to  the  energy  levels  of  each  element  in  the 
specimen  being  irradiated.  The  relative  intensities  of  the  various  peaks  will 
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thus  be  a  measure  of  the  chemical  composition  of  the  specimen. 

The  binding  energy  is  actually  a  measure  of  the  potential  around  the 
atom  of  interest.  The  same  atom  in  different  chemical  environments  will  have 
a  different  potential  depending  on  the  valence  charge  distribution.  Binding  to 
more  electronegative  atoms  will  result  in  charge  transfer  to  the  more 
electronegative  atoms  from  the  atom  of  interest,  decreasing  the  valence 
charge  density  and  thus  increasing  the  observed  binding  energy.  The 
binding  energies  of  the  ohotoeiectron  peaks  are  thus  not  only  a  measure  of 
the  elemental  composition  of  the  specimen,  but  also  are  a  measure  of  the 
oxidation  states  of  each  of  the  constituent  elements.  Differences  in  binding 
energies  which  are  due  to  differences  in  the  chemical  environment  are 
termed  chemical  shifts.  Changes  in  the  valence  charge  distribution  can  also 
arise  from  structural  differences,  which  may  also  result  in  observable  binding 
energy  shifts.  One  example  of  such  structure  induced  charge  transfer  (SICT) 
occurs1'3  in  the  previously  mentioned  case  of  SiC^/Si,  in  which  the  structure 
of  the  interfacial  SiC>2  differs  from  that  of  the  bulk  SiC>2. 

It  is  to  be  noted  that  XPS  is  a  surface  sensitive  technique,  since  the 
photoelectron  inelastic  mean  free  path  (the  path  length  over  which  the 
fraction  of  photoeiectrons  which  are  not  inelastically  scattered  is  1/e)  is 
typically  in  the  range  5-50  A,  depending  upon  the  material  and  the  kinetic 
energy  of  the  emitted  electron.  This  surface  sensitivity  can  be  exploited  to 
obtain  compositional  depth  profiles.  This  can  be  accomplished  by  varying  the 
photon  energy,29  and  hence  the  photoelectron  kinetic  energy  and. the  mean 
free  path.  Alternatively,  the  electron  take-off  angle  can  be  varied.30'34  In  this 
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work,  sequential  chemical  etching  has  been  used  to  exoose  the  insulator  film 
at  different  depths  below  the  initial  surface  to  the  vacuum.  Details  of  the 
etching  technique  can  be  found  elsewnere.’ -35  The  intensity  analysis 
necessary  to  obtain  a  compositional  depth  profile  using  :n:s  approach  was 
developed17  as  a  pad  of  this  task,  and  is  a  generalization  of  a  previously 
published  technique.36  The  intensity  analysis  is  discussed  in  detail  in 
Appendix  A. 

In  some  cases,  it  is  desirable  to  study  the  insulator/Si  interface  from  the 
Si  side.  Such  a  situation  occurs  in  device  structures  whicn  have  been 
metallized,  since  the  metal  must  be  removed  uniformly  without  affecting  the 
insulator  fiim,  which  is  difficult.  Another  such  situation  occurs  in  the  study  of 
the  buried  insulator  layer  in  SOI  structures.  Finally,  verification  of  the 
interfacial  composition  may  be  required  if  it  is  suspected  that  the  chemical 
etch  resulted  in  preferential  etching,  yielding  results  which  are  artifacts  of  the 
etching  technique.  In  all  of  these  cases,  an  etch  of  Si  which  does  not  affect 
the  insulator  film  is  required.  XeF2  gas-phase  etching  of  Si  has  been 
developed37-38  as  part  of  this  task,  and  has  been  used  in  parts  of  this  work. 

While  XPS  has  proven  to  b€  a  powerful  technique  for  determining 
chemicai  compositions,  it  is  not  the  best  technique  for  all  circumstances. 
Chemical  shifts  can  be  small,  signals  may  be  weak,  and  it  is  difficult  to 
unambiguously  determine  structural  differences  from  XPS  spectra  alone.  In 
addition,  the  surface  sensitivity  of  XPS  can  be  a  disadvantage  if  a  bulk 
property  of  a  material  is  of  interest,  or  if  an  interface  to  be  studied- is  beyond 
the  detection  depth  of  XPS  and  there  is  no  suitable  chemical  etch  to  expose 
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the  interface.  In  such  circumstances,  other  techniques  must  be  used. 

In  this  work,  spectroscopic  (wavelength  dispersive)  ellipsometry  has 
been  used  :n  conjunction  with  XPS.  Ellipsometry  involves  measurement  of 
the  polarization  state  of  iight  reflected  from,  a  specimen  when  the  ooiarization 
state  of  tne  incident  light  is  known.  Reviews  of  the  technique  and  its  concepts 
and  applications  are  available  elsewhere,39'42  so  only  a  brief  description  of 
the  concepts  is  presented  below.  The  spectroscopic  ellipsometer  used  in  this 
work  is  similar  to  that  described  by  Aspnes  and  Studna.43  and  is  descnbed 
elsewhere.44 

With  the  incident  and  reflected  beams  defining  the  plane  of  incidence, 
the  polarization  state  of  the  light  can  be  determined  in  tornr  of  a  component 
linearly  polarized  in  the  plane  of  incidence  (p- polarized),  and  a  component 
linearly  polarized  perpendicular  to  the  plane  of  incidence,  in  the  plane  of  the 
sample  surface  (s-polarized,  from  the  German  senkrecht  for  perpendicular). 
The  ratio  of  the  electric  vectors  of  the  reflected  and  incident  beams  is  defined 
as  the  reflection  coefficient  R.  The  complex  ratio  of  the  reflected  p-  and 
s-polarized  light  components  is  related  to  the  measured  ellipsometric 
parameters  tan  y  and  cos  A  via 

Rp(E)/Rs(E)  =  p(E)  =  tan  y(E)  exp[iA(E)],  0  -1 ! 

where  E  is  the  photon  energy.  The  material  property  to  which  p(E)  is  related 
is  the  dielectric  function  e(E),  which  is  defined  as  the  permittivity  relative  to  the 
free  space  permittivity  a0. 
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The  region  of  the  sample  which  contributes  to  the  ellipsometnc 
measurement  depends  upon  how  absorbing  the  sample  is  at  the  plioion 
energy  used  in  the  measurement.  For  Si  and  GaAs,  the  absorption  depth  (the 
distance  over  which  the  light  intensity  is  aitenuatec  by  a  factor  cx  e 
analogous  to  the  photoelectron  mean  free  path  in  XPSf  in  the  visible  - 
ultraviolet  energy  range  can  De  as  low  as  50  -  100  A,  as  large  as  thousancs 
of  A  near  the  bandgap,  and  much  larger  below  the  bandgap  where  the 
materials  are  transparent.  In  its  ability  to  "see"  through  thick  fiims  in  an 
appropriate  photon  energy  range,  SE  is  thus  complementary  to  XPS. 
However,  SE  can  also  be  surface  sensitive,  since  reflectivity  can  be 
sinnific.antly  affected  by  the  surface  condition  of  the  sample. 

Spectroscopic  ellipsometry  differs  from  conventional  null  ellipsometry  in 
its  use  of  variable  wavelength  light,  and  its  use  of  a  continuously  rotating 
polarizer  or  analyzer  together  with  Fourier  analysis  of  the  time-dependent 
light  intensity  at  the  detector,39'44  making  a  compensator  unnecessary.  The 
Fourier  analysis  yields  the  ellipsometric  parameters  tan  y  and  cos  A.  The 
pseudodielectric  function  (the  measured  dielectric  function  assuming  a 
mathematically  abrupt  solid-ambient  interface),  <e(E)>,  can  be  calculated 
from  p(E)  if  the  sample  can  be  considered  as  a  semi-infinite  medium  of  a 
single  material,  rather  than  a  multilayer  structure  of  more  than  one  material. 
The  semi-infinite  medium  may  itself  consist  of  a  physical  mixture  of  more  than 
one  material,  in  which  case  the  mixture  can  be  modeled  by  using  an  effective 
medium  approximation45  46  to  calculate  the  effective  dielectric  function  of  the 
mixture.  The  composition  of  the  mixture  is  varied  to  achieve  good  agreement 
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with  the  data.  The  fraction  of  each  material  in  the  mixture  can  thus  oe 
obtained  numerically.  The  density  of  a  material  is  a  bulk  physical  property 
which  can  be  obtained  by  modeling  the  material  as  the  same  material  ;r. 
some  known  standard  state  with  an  admixture  of  voids. 

Spectroscopic  eliipsometric  data  obtained  from  multilayer  structures  is 
analyzed  in  a  manner  similar  to  that  discussed  above,  though  the  relevant 
equations  are  complicated  by  interference  between  light  reflected  from  the 
different  interfaces  in  the  multilayer  structure,  and  by  multiple  internal 
reflections  within  each  layer.  Physical  properties  which  can  be  obtained 
using  multilayer  analysis  include  the  degree  of  surface  roughness  (by 
modeling  the  surface  as  a  layer  consisting  of  a  physical  mixture  of  the 
underlying  material  and  voids),  interfacial  roughness  (by  modeling  the 
interface  as  a  layer  consisting  of  a  physical  mixture  of  the  materials  on  either 
side  of  the  interface),  and  the  film  thickness  and  composition. 

In  the  discussion  above,  the  use  of  SE  to  determine  physical  properties 
of  materials,  i.  e.  properties  which  would  be  considered  as  structure-related  in 
the  context  of  the  earlier  discussion,  has  been  briefly  considered.  Chemical 
information  can  also  be  ODtained  using  SE.  For  example,  the  chemical 
composition  of  the  SiC^/Si  interface  has  been  measured47  using  calculated 
values48  of  the  dielectric  function  of  SiOx.  However,  SE  has  most  commonly 
been  used,  and  is  best  suited,  for  the  determination  of  structure-related 
properties  of  materials,  and  is  in  this  respect  complementary  to  XPS. 

The  remainder  of  this  report  is  organized  as  follows.  Studies  of  metal- 
nitride-oxide-semiconductor  structures  are  described  in  Chapter  2,  studies  of 
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SiOxNy/Si  formed  by  the  thermal  nitndation  ot  S 1 0 2 ' S i  are  discussed  in 
Chapter  3,  and  studies  of  semiconductor-on-insulator  structures  are 
described  in  Chapter  4.  The  conclusions  are  contained  in  Chapter  5, 
followed  by  the  bibliography  ana  appendix. 
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Chapter  2 

Studies  of  Metal-Nitride-Oxide-Semiconductor  (MNOS) 

Structures 


As  previously  mentioned,  the  study  of  CVD  S13N4  and  its  interface  with 
S1O9  m  MNOS  structures  was  the  primary  work  of  the  previous  task,  and  the 
results  are  described  elsewhere.4'9  Some  foliow-on  work  has  been  done  in 
this  task,  concerning  the  effects  of  H2  annealing  on  the  S13N4  anc  the 
S'3N4  S1O2  interface,  and  the  results  are  described  herein. 

Two  samples  were  obtained  from  Westmghouse  which  consistec  o* 
aporoximately  450  A  of  S13N4  deposited  on  Si  wafers  whicn  had  native 
oxides.  One  wafer  was  annealed  in  H2,  while  the  other  was  used  as  a  control 
and  was  not  annealed.  Both  wafers  were  then  metallized  with  Al  to  form  the 
MNOS  structure  and  shipped  to  JPL.  The  metallized  surface  was  bonded  to  a 
Au  platen  with  In  and  the  Si  wafer  was  etched  away  with  XeF2  gas.  a 
technique  described  elsewhere. 3'3S  Briefly,  XeF2  reacts  spontaneously  witn 
Si  to  form  Xe  and  SiF4>  but  does  not  react  with  S1O2  or  S13N4  unless 
activated  with  a  beam  of  electrons,  ions,  or  UV  photons.  XeF2  is  a  soiid 
which  has  a  room  temperature  vapor  pressure  of  4.5  torr.  Since  the  reaction 
products  of  XeF2  with  Si  are  gases,  a  simple  vacuum  system  with  a  valve  to 
separate  the  XeF2  holding  chamber  and  the  sample  chamber  can  be  used  to 
accomplish  the  etcning.  Figure  2.1  shows  a  schematic  diagram  of  the  XeF2 
etching  system  used  in  this  work.  Using  this  technique,  the  native  S1O2  Si 
'nterface  was  exposed.  XPS  spectra  were  recorded,  then  both  samples  were 
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VAlVt 


e  system  used  for  gas  phase  XeF- 


etcnea  with  1:10  HF:E*OH  to  expose  the  Si02'Si3N4  interface  arc  X 
soectra  were  again  accumulated. 

Prior  to  the  HF:EtOH  etch,  the  Si  2d  soect-a  from  ooth  samples  srcv 
two  distinct  Dea*s  separatee  by  .5  eV.  t re  same  separatmm  as  : 
coserved'3  between  SIO2  and  S13N4.  In  addition,  the  erergy  seoa-a 
between  the  Si  2 p  peak  at  higher  energy  and  the  O  is  oea*  s  -29  ~ 
wmch  is  within  experimental  error  of  the  ^-23.6  eV  expected*  fcr  S : 0 ~ 
confirms  that  some  SiO?  remains  intact  at  the  S102-S.  inter'ace  ah~r 
S:3N4  cebositicn,  as  previously  reported'-9  from  studies  c*  sa~c  es 
"30  been  ceoth  profiled  from  the  front  sice.  The  earner  firemgs' are  t 
not  due  to  possible  artifacts  of  the  etching  technique.  In  this  wor-:.  :  - 
‘ound  that  the  two  samples  had  native  oxices  wmch  were  semewa:  c  "e' 
'■n  thickness,  with  the  H2 -annealed  sample  having  a  t.mnner  -at  ve  ox  oe 
A  vs.  20  A  fcr  the  control  sample.  These  thiCKnesses  were  est  ~a-ec 
assuming  a  30  A  ohotoelectron  mean  free  path  in  both  5  C?  arc  S  -;N4  a 
using  intensity  analysis  as  described  :n  Appendix  A  arc  eisewnere  *' 

After  the  HF'EfOH  etch,  the  Si  2 p  spectra  from  the  two  sauc  es  ?-ov 
substantial  differences.  The  P^-anneaied  sample  again  snowed  two  ce: 
with  the  energy  separation  reduced  to  1.2  eV.  indicative  of  an  cxy-mo- 
Only  one  Si  2 p  peak  was  apparent  for  the  control  sample,  c~: 
analysis49  suggested  the  presence  of  twe  peaKS  wtn  a  prcacer  mesm 
and  wider  energy  separation  than  was  observed  for  the  H~-ar','ea  ec  sa~ 

The  data  can  be  explained  in  the  ‘allowing  manner  T~e  H-  ‘or”s  5  - 
groups  m  the  oxide  and  Si-NH?  m  the  nitride  ;n  a  manner  si”'  ar  *o  * 


proposed'8  for  the  thermal  nitndation  of  SiC>2-  The  hydroxyl  and  amine 
groups  are  mobile  at  high  temperature,  resulting  in  interdiffusion  at  the 
SiOj'S^N^  interface  with  subsequent  oxvnitride  formation.  This  explains  the 
iower  energy  seoarat.on  between  the  two  Si  2 p  peaKS  after  ho  anneaiing 
The  narrower  lineshape  observed  after  H2  anneali  g  may  be  a  result  of  strain 
relaxation.  As  previously  reported,''3  strain  near  the  S1O2  Si  interface  results 
in  a  different  distribution  of  Si-O-Si  bond  angles  than  that  which  occurs  in 
bulk  SiOo,  with  an  accompanying  charge  transfer  which  shifts  the  binding 
energy  of  the  Si  2 p  peak.  Since  different  bond  angles  result  in  different 
binding  energies,  the  effect  is  an  apparent  broadening  in  the  Si  2 p  lineshape. 
The  breaking  of  strained  Si-O-Si  bonds  by  the  H2  to  form  Si-OH  would  be 
exoected  to  relieve  the  strain  and  allow  the  unbroken  Si-O-Si  to  relax  closer 
to  the  144”  bond  angle  which  is  predominant  in  bulk  SiO £.  The  reduction  in 
the  number  of  Si-O-Si  bond  angles  with  higher  and  lower  values  than  144” 
would  thus  be  expected  to  narrow  the  Si  2 p  lineshape. 

In  summary,  the  effect  of  H2  annealing  appears  to  be  to  make  the 
Si3N4/Si02  interface  less  abrupt  by  causing  Interdiffusion  of  silanol  and 
silamine  groups  and  subsequent  oxynitride  formation,  and  also  to  relieve 
strain  at  the  Si02//Si  interface. 
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Chapter  3 

Studies  of  SiOxNy/Si  Formed  by  the  Thermal 
Nitridation  of  Si'02/Si 

I.  Kinetics  of  Nitridation:  Background 

S1O2  films  on  Si  have  been  extensively  studied1-3  and  are  known  to  be 
chemically  homogeneous.  The  reaction  of  Si02  with  NH3  results  in  the 
formation  of  an  oxynitride.  This  material  may  have  chemical 
inhomogeneities,  and  the  fact  that  the  progenies  of  the  initial  S1O2  film  are 
well  understood  may  make  it  possible  to  gain  a  detailed  understanding  of  the 
formation  of  these  inhomogeneities  by  studying  and  modeling  the  kinetics  of 
nitridation. 

Silicon  oxynitride  (SiOxNy),  formed  by  the  thermal  nitridation  in  NH3  of 
Si03  films  on  Si,  is  a  promising  dielectric50'58  for  such  applications  as  very 
large  scale  integration  (VLSI),  floating  gate  electrically  erasable 
programmable  read-only  memory  (EEPROM)  devices,  and  radiation  hard 
devices.  The  thin  gate  or  tunnel  insulators  required  for  VLSI  and  floating  gate 
EEPROM  aevices  can  be  a  source  of  serious  problems  if  the  insulator  is  the 
conventional  thermally  grown  Si02-  These  problems  include50"58  low 
resistance  to  dopant  and  impurity  diffusion,  high  electric  field  instabilities,  hot 
electron  effects,  radiation  effects,  and  slow  trapping.  SiOxNy  is  in  many  ways 
superior  to  Si03  in  the  applications  mentioned  above.  Specifically,  it  has 
been  shown50'56  that  thermally  nitrided  Si03  has  an  increased  resistance  to 
high  field  instability,  forms  a  better  diffusion  barrier,  is  less  reactive  to  gate 
electrodes,  is  oxidation  resistant,  reduces  degradation  due  to  carrier  injection. 
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and  nearly  eliminates  formation  of  ionizing  radiation-induced  interface  states. 
S13N4  exhibits  many  of  these  same  advantages.  Direct  thermal  mtridation  of 
Si  has  been  shown  to  be  possible,  but  the  thickness  of  the  film  is 

self-limited59  at  about  30  A.  In  addition,  deposited  S13N4  has  been  utilized  in 
MNOS  (metai-nitnde-oxide-semiconductcr)  structures.  However,  this  duai 
dielectric  structure  exhibits53  high  trapping  efficiency  in  the  nitride  and  at  the 
nitride/oxide  interface,  making  it  unsuitable  for  some  applications.  Thermally 
nitrided  SiOg  has  the  advantages  of  not  having  the  dielectric/dielectric 
interface  of  MNOS  structures,  and  the  thickness  is  determined  by  the  initial 
oxide  thickness. 

The  potential  use  of  SiOxNy  in  technologically  important  device 
apDlications  has  prompted  several  studies10'18'50'53’60'69  of  the  chemical 
structure  of  S1O2  films  which  have  been  thermally  nitrided  in  NH3.  Early 
studies50'53'60*65  were  done  for  varying  nitridation  conditions,  which  (not 
surprisingly)  resulted  in  reported  N  distributions  which  differed  from  study  to 
study.  There  was  general  agreement  that  the  nitrogen  concentration  in  the 
films  is  higher  at  the  surface50'52'63'65  and  near  the  SiOxNy/Si 
interface62'54'65  than  in  the  bulk  of  the  oxynitride  film,  although  two  studies50'03 
found  direct  evidence  only  for  the  higher  surface  nitrogen  concentration,  and 
one  study60  found  the  nitrogen  distributed  uniformly 
throughout  the  film.  Recently,  it  was  reported13  that  the  peak  in  the  interfacial 
nitrogen  concentration  can  occur  away  from  the  SiOxNy/Si  interface  itself. 

Some  of  the  differences  noted  above  may  be  explainable  by  the  different 
sensitivities  of  the  techniques  employed,  since  the  AES  sputter  depth  profiling 
•ypically  employed  is  known  to  broaden  the  interface,  and  thus  resolution  of 
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possible  differences  in  N  distributions  in  the  critical  interfacial  region  is 
difficult.  However,  some  of  the  differences  are  almost  surely  related  to  the 
kinetics  of  the  mtridation  process  and  differences  in  nitridation  conditions. 
The  findings  discussed  above  have  prompted  several  recent  studies14'18-67'63 
of  the  kinetics  of  the  nitridation  process.  The  work10-18  described  in  detail  in 
this  chapter,  together  with  work67'69  done  independently  and  in  parallel  at 
Stanford,  represents  the  first  comprehensive  attempt  to  determine  the 
dependence  of  the  N  distribution  on  nitridation  time  and  temperature.  In 
addition,  this  work10'18  is  the  first  to  use  a  technique,  high  resolution  XPS  in 
conjunction  with  relatively  benign  chemical  etching,  which  is  capable  of  a 
detailed  study  of  the  interfacial  nitrogen  distribution.  The  reliable 
determination  of  the  interfacial  nitrogen  distributions  has  allowed,  for  the  first 
time,  a  detailed  picture  of  the  mechanism  and  kinetics  of  the  nitridation 
process  to  be  proposed  and  computer  simulations18  of  the  nitridation  kinetics 
carried  out  as  part  of  this  work.  The  kinetic  modeling  also  shows  that  the  role 
of  interfacial  strain,  discussed  in  detail  below  and  in  section  IV,  must  be  taken 
into  account  to  explain  the  observed  interfacial  nitrogen  distnbutions. 

The  dependence  of  electrical  characteristics  on  the  nitndation  conditions 
has  also  been  reported  in  the  literature,  as  described  in  detail  below  and  in 
section  III.B.  Thus,  as  a  further  check  of  our  proposed  model  for  the  kinetics  of 
nitridation  we  have  examined  these  data  and  shown,  for  the  first  time,  a 
correlation  with  the  nitrogen  distributions  observed  in  this  work.  This 
correlation  is  explained  within  a  model15  in  which  the  defect  formation 
kinetics  are  postulated  to  be  strain-dependent,  and  correlated  to  the  strain 
dependent  kinetics  of  the  nitridation  process  itself.  This  work  thus  presents 
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the  first  comprehensive  explanation  of  the  interdependence  of  the  processing 
conditions,  electrical  characteristics,  and  chemical  structure  of  nitrided  oxides 
within  a  framework  which  takes  into  account  known  properties  of  the  initial 
S1O2  film. 

Of  the  electrical  studies  alluded  to  in  the  preceeding  paragraph,  one 
study52  reported  that  the  interface  fixed  charge  density  correlates  with  the 
amount  of  nitrogen  at  the  interface,  while  another  reported53  that  the 
magnitude  of  the  flatband  voltage  shift  |AV FB{  --  typically  interpreted  as  a 
measure  of  the  fixed  charge  in  the  dielectric  --  is  reduced  by  the  presence  of 
nitrogen  at  the  interface.  A  more  recent  report57  showed  that  for  a  nitridation 
temperature  of  900‘C,  | A V FB |  increases  monotonically  with  increasing 
nitridation  time,  but  for  a  nitridation  temperature  of  950*C,  |A V FB |  initially 
increases,  then  decreases  with  increasing  nitridation  time.57  At  nitridation 
temperatures  >  1000'C,  |AVFB|  decreases  monotonically57  for  nitridation  times 
beyond  30  minutes.  Similar  observations  on  the  behavior  of  |AVFB|  have 
been  made  by  other  researchers.70 

Chen  et  at.  suggested57  that  the  lower  values  of  |AVFB|  for  SiOxNy/Si 
formed  at  high  nitridation  temperatures  (>  1000'C)  may  be  due  to  a  small 
solubility  of  hydrogen  --  which  was  postulated  to  be  involved  in  the  formation 
of  defects  --  at  elevated  temperatures.  However,  the  solubility  of  H2  in  fused 
silica  is  nearly  independent  of  temperature.71  The  hydrogen  may  also  be  in 
the  form  of  H2O,  which  reacts  with  SiC>2  to  form  Si-OH  bonds,  but  the 
solubility  of  H2O,  and  hence  the  concentration  of  Si-OH  bonds,  in  fused  silica 
is  also  not  a  strong  function  of  temperature.72-73  One  may  thus  conjecture  that 
the  same  holds  true  for  the  solubility  of  NH3  in  Si02-  It  is  more  likely  that  the 
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formation  of  the  defects  giving  rise  to  |AVFBj  is  related  to  the  nitridation 
kinetics.  In  this  work,15  a  correlation  between  the  nitrogen  distribution  and  the 
behavior  of  |AVFB;  will  be  identified  and  an  explanation  will  be  proposed. 

In  order  to  understand  the  underlying  reasons  for  the  measured  nitrogen 
distributions,  we  have  developed18  a  kinetic  model  which  simulates  the 
nitridation  process.  In  any  attempt  to  model  the  nitridation  kinetics,  the  initial 
structure  of  the  Si02  film  must  first  be  understood.  In  particular,  the  possibility 
that  strain  in  the  Si02  film  affects  the  nitridation  kineses  must  be  considered. 
It  is  known  that,  upon  oxidation  of  Si  to  form  SiC^.  a  large  molar  volume 
change  oroduces  compressive  stress  in  the  oxide  in  the  plane  of  the  Si 
surface,  which  relaxes  by  viscous  flow  in  the  growth  direction.74'75  At 
oxidation  temperatures  below  950‘C,  the  relaxation  rate  is  slower  than  the 
oxidation  rate,74  resulting  in  a  buildup  of  intrinsic  stress  in  the  oxide  and 
denser  oxides  at  low  oxidation  temperatures.  This  temperature  will  be 
referred  to  as  the  viscous  flow  temperature  Tv.  The  stress  and  densification  in 
the  oxide  at  low  growth  temperatures  may  affect  the  nitridation  kinetics.  This 
first  type  of  strain  is  macroscopic  in  nature,  affecting  the  entire  oxide  film. 

Another  source  of  strain  is  the  lattice  mismatch  between  the  Si  and  the 
Si02,  and  is  localized  w'thin  30  A  of  the  Si02/Si  interface.1'3  This  interfacial 
strain  is  present  even  in  high  temperature  oxides1'3’36'76'77  which  do  not 
exhibit  macroscopic  stress.74  This  strained  interfacial  region  is  characterized 
by  a  higher  proportion  of  4-membered  rings  (4  Si  tetrahedra  with  bridging 
oxygens  in  a  ring  structure)  compared  to  the  unstrained  bulk  SiC^.  which 
consists  primarily  of  6-membered  nngs.  The  effects  of  the  strain  would  not  be 
expected  to  be  confined  to  the  SiC>2  side  of  the  interface,  but  would  extend  a 
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few  A  into  the  Si  substrate  as  well.  Interfacial  strain  has  been  invoked  to 
explain'8  the  observed  U-shaped  density  of  states  in  tne  Si  band  gap  at  the 
S i / S i O 2  interface,  and  in  general  in  lattice  mismatched  systems,  and  to 
explain  the  generation76-'9'80  and  diffusion77  of  defecrs  in  electron-  irradiated 
SiS>2  Si.  It  will  be  shown  that  the  measured  nifogen  attributions  are 
consistent  with  the  kinetic  model  providec  the  influence  of  the  interracial 
strain  on  the  nitridation  kinetics  is  explicitly  accounted  for. 

The  overall  objective  of  this  work10'18  has  been  to  determine  the  kinetics 
and  mechanism  of  the  nitridation  process.  These  studies  have  utilized  high 
resolution  XPS  in  conjunction  with  nondestructive  chemical  depth 
profiling13-14-16-17  to  measure  nitrogen  distributions,  kinetic  simulations  of  the 
nitridation  process14-16*18  to  understand  the  factors  giving  rise  to  the 
measured  nitrogen  distributions,  and  spectroscopic  ellipsometry  for  possible 
complementary  information.  XPS  is  not  well  suited  for  determining  physical 
properties  such  as  void  fractions,  interfacial  roughness,  or  lateral 
nonuniformities.  Spectroscopic  ellipsometry  can,  in  principle,  provide 
information  on  these  physical  properties.  Together,  these  techniques  give 
insight  into  how  processing  conditions  affect  the  chemical  composition  of  the 
nitnded  oxides,  and  provide  a  basic  understanding  of  the  nitridation  process 
and  the  role  of  strain  in  the  nitridation  kinetics.  In  addition  to  determinina 
chemical  compositions,  defect  distributions  have  also  been  deduced,  both  by 
correlating’5  chemical  compositions  with  electrical  characteristics,  and 
through  use16-17  of  a  fluorine  marker,  which  is  postulated  to  bind  to  defect 
sites.  Some  understanding  of  the  relationship  between  processing 
conditions,  chemical  composition,  and  electrical  characteristics  can  thus  be 
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The  experimental  conditions  will  oe  described  in  Section  II.  The  XPS 
measurements  will  be  described  in  Section  III,  with  the  measured  nitrogen 
distributions  in  mtrided  oxides  presented  in  Section  III. A  and  the  defect 
distributions  deduced  from  electrical  and  chemical  data  presentee  in  Section 

III. B,  together  with  a  model  to  explain  the  electrical  data.  Section  IV  contains 
the  kinetic  simulations,  with  the  nitridation  process  described  in  Section  IV. A, 
the  kinetic  model  presented  in  Section  IV. B,  the  choice  of  kinetic  parameters 
described  in  Section  IV. C,  and  the  simulation  results  presented  in  Section 

IV.  D.  The  spectroscopic  ellipsometry  results  will  be  described  in  Section  V, 
and  Section  VI  contains  the  summary  and  conclusions. 

II.  Experimental 

Si  (100)  wafers  were  thermally  oxidized  in  dry  O2  to  a  thickness  of 
approximately  100  A,  then  heated  in  pure  NH3  at  atmospheric  pressure. 
Oxidation  temperatures,  nitridation  temperatures,  and  nitridation  times  are 
shown  in  Table  I.  The  XPS  spectrometer  used  in  this  work  is  a  modified 
HP5950A1  utilizing  monochromatized  AIKa  x-rays  (1486.6  eV).  The  sample 
introduction  area  of  the  spectrometer  is  enclosed  by  a  dry  box  in  which  an  N3 
ambient  is  maintained  in  order  to  prevent  exposure  of  the  sample  to 
atmospheric  oxygen  and  water  vapor  after  it  has  been  etched.  The  samples 
were  sequentially  chemically  etched  in  the  dry  box  with  1:10  HF:EtOH 
dispensed  with  a  microliter  pipette  onto  a  sample  spinning  at  3600  rpm.  The 
EtOH  used  was  U.S.P.  punctilious  grade  ethanol.  Details  of  the  etching 
technique  and  its  advantages  have  been  discussed  elsewhere.1-35  This 


Table  I.  Summary  of  the  oxidation  and  nitridation  conditions  for  the 
samples  used  in  this  work. 


3mple  # 

Initial  Oxide 
Thickness 

(A) 

Oxidation 

Temperature 

(°C) 

Nitridation 

Temperature 

(°C) 

Nitridation 

Time 

(min) 

1 

100 

1000 

1000 

10 

2 

100 

1000 

1000 

30 

3 

100 

1000 

1000 

1 20 

4 

100 

825 

1000 

240 

5 

100 

1000 

1150 

5 

6 

100 

1000 

1150 

10 

7 

100 

825 

800 

60 

8 

100 

1000 

800 

60 

9 

100 

1000 

800 

240 

10 

1030 

1000 

800 

60 

1 1 

1030 

1000 

1000 

60 
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chemical  depth  profiling  technique  preserves  chemical  and  structural 
information  which  would  be  lost1  if  ion  sputtering  techniques  were  employed. 

The  base  pressure  of  the  analysis  chamber  was  2-3  x  10‘10  Torr  during 
accumulation  of  SDectra.  After  each  etcn  step.  XPS  spectra  were  recorded 
and  the  Si  2 p,  0  is,  N  Is,  and  F  Is  photoelectron  peak  intensities  were 
measured  by  fitting  Voigt  profiles  (convolution  of  a  Lorentzian  function 
representing  core  hole  lifetime  broadening  and  a  Gaussian  function 
representing  instrumental  and  other  broadening  mechanisms)49  to  the  peaks 
and  correcting  the  intensities  using  Scofield's  calculated  photoionization 
cross  sections.31  These  intensities  were  used  to  calculate  the  relative 
nitrogen  concentration  [N]/([0]+[N])  by  the  technique  described  in  Appendix  A. 
Using  this  technique,  the  actual  film  composition  as  a  function  of  depth  is 
obtained,  unbiased  by  the  averaging  effect  of  the  relatively  long 
photoelectron  mean  free  path. 

III.  XPS  Measurements 

The  nitrogen  distribution  measured  in  nitrided  oxides  will  be  described  in 
section  III. A.  In  section  II  1.3 ,  the  defect  distributions  deduced  from  electrical 
and  chemical  data  are  presented  and  correlated  to  the  observed  nitrogen 
distributions.  A  model  is  proposed  to  explain  the  correlation  between  the 
electrical  data  and  nitrogen  distributions. 

III.  A.  Nitrogen  Distributions 

Figure  3.1  shows  the  N  1  s  and  O  Is  photoelectron  peak  intensities,  lN1s 
and  IQ1  s,  for  sample  1 ,  which  has  been  nitrided  at  1 000’C  for  1 0  minutes.  The 
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intensities  have  been  corrected  for  the  different  photoionization  cross 
sections.8'  The  inhcmogeneous  nitrogen  distribution  in  the  mtriaed  oxide  is 
clearly  evident,  with  the  nitrogen  signal  high  relative  to  the  oxygen  signal  in 
the  surface  and  intedacial  regions.  This  is  consistent  with  the  earliest 
studies50-52  53'65  of  the  distribution  of  r-trogen  in  nitrided  oxides.  With 
increased  nitridation  time,  however,  lN,s  decreases  faster  than  l0,s  in  the 
region  closest  to  the  interface.  This  is  consistent  with  a  later  report13  that  the 
peak  interfacial  nitrogen  concentration  occurs  approximately  20  A  away  from 
the  interface.  These  trends  are  more  apparent  when  the  ratio  lN1s/(l0is+  lN.s; 
is  plotted,  as  shown  in  Fig.  3.2.  However,  this  is  not  a  measure  of  the  actual 
nitrogen  distribution  m  the  film,  since  a  correction  for  the  difference  in  the 
mean  free  pathn  ot  N  Is  and  O  Is  photoelectrons  (due  to  the  different  kinetic 
energies)  has  not  been  included.  In  addition,  the  averaging  effect  of  the 
relatively  long  mean  free  paths  has  not  been  considered.  These  corrections 
have  been  made  using  the  technique  described  in  Appendix  A  for  the  data 
shown  in  Fig.  3.3,  which  shows  the  nitrogen  distributions  for  samples  1,  2. 
and  3.  For  all  three  samples,  the  oxidation  and  nitridation  temperatures  are 
both  1000°C,  and  the  nitridation  times  range  from  10  to  120  minutes. 

For  a  nitridation  temperature  of  1000°C,  the  nitridation  time  dependence 
of  the  nitrogen  distribution  can  be  summarized  as  follows.  Nitrogen  initially 
incorporates  rapidly  at  the  interface  until  a  saturation  concentration  is 
reached,  as  seen  in  Fig.  3.3(a).  As  the  nitridation  proceeds,  the  surface 
nitrogen  concentration  also  saturates,  nitrogen  incorporation  continues  in  the 
bulk  of  the  film,  and  the  interfaciai  nitrogen  peak  occurs  away  from  the 
interface,  as  seen  in  Fig.  3.3(b).  With  increasing  nitridation  time,  nitrogen 
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Figure  3.3.  Distribution  of  N  in  SiOxN  produced  from  a  100  A-thick  SiCs  S'm 
grown  m  dry  02  at  1000'C,  then  mtrided  at  1000’C  for  fa)  10  minutes 
(b)  30  minutes,  (c)  120  minutes. 
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Distance  from  SiO.N  /Si  Inter  luce  (A) 


incorporation  continues  in  the  bulk  of  the  dielectric,  as  seen  in  Fig.  3.3(c,'  and 
in  the  data  from  sample  4,  which  has  been  previously  reported.'3 
Qualitatively  similar  results  have  been  reported67'69  for  400-A-thick  oxides 
mtrioed  at  110C3C.  Quantitative  differences  can  be  attributed  to  the  different 
initial  oxioe  thickness  anc  higner  nitridation  temperature  used  in  that 
work.0,  '°9 

F  ]ure  3.4  shows  the  nitrogen  distribution  in  samples  5  and  6,  which 
have  been  mtrided  at  1150°C.  Figures  3.4(a)  and  3.4(b)  are  similar  to 
Figures  3.3(b)  and  3.3(c),  except  the  peak  interfacial  nitrogen  concentration 
has  saturated  at  a  lower  value,  although  it  occurs  essentially  the  same 
distance  from  the  interface.  The  ratio  [N]/([0]  +  [N])  is  seen  to  have  a 
maximum  value  near  0.55  in  the  interfacial  region  for  a  nitridation 
temperature  of  1 000'C,  and  near  0.45  for  a  nitridation  temperature  of  1 1 50'C. 
The  decreased  nitrogen  incorporation  in  the  interfacial  region  at  elevated 
temperatures  may  partially  explain  the  failure  to  observe50  an  interfacia! 
nitrogen  buildup  for  a  500  A  thick  oxide  nitrided  at  1200  C. 

Figure  3.5  shows  the  nitrogen  distributions  in  samples  7,  8,  and  9,  which 
have  been  mtrided  at  800°C.  Samples  7  and  8  differ  only  in  the  oxidation 
temperature.  The  purpose  of  this  comparison  is  to  determine  the  effect  of 
macroscopic  intrinsic  stress  on  the  nitridation  kinetics.  The  one  hour 
nitridation  time  for  samples  7  and  8  is  insufficient74  for  the  macroscopic  stress 
to  relax  by  viscous  flow  at  the  nitridation  temperature  of  800°C,  so  the  low 
temperature  oxide  would  retain  the  stress  acquired  during  oxidation.  It  can 
be  seen  in  Fig.  3.5(a),  however,  that  the  nitridation  has  not  been  affected  by 
the  oxidation  temperature.  This  may  be  because  the  effect  of  macroscopic 
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Figure  3.4.  Distribution  of  N  in  SiO  N  produced  by  nitriding  a  100  A-thick 
Si02  film  at  1 1 50'C  for  (a)  5  minutes,  and  (b)  1 0  minutes. 
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Figure  3.5.  (a)  Comparison  of  the  distribution  of  N  in  SiO  N  films  produced 

from  100  A-thick  Si02  films  grown  in  dry  02  at  825’Cf  (solid  circles)  and 
1000’C  (open  circles),  then  nitrided  at  800'C  for  60  minutes,  (b) 
Distribution  of  N  in  a  SiO  N  film  produced  from  a  100  A-thick  SiO,  film 
nitrided  at  800'C  for  240  mihutes.  2 
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stress  on  the  mtridation  is  too  small  to  measure,  or  because  the  low 
temperature  oxide  has  relaxed  faster  than  expected  due  to  the  formation  of 
Si-OH  and  Si-NHX  in  the  oxynitride.  It  has  previously  been  reported74  that 
relaxation  is  faster  in  a  h^O-containing  ambient,  which  reacts  with  SiO?  to 
form  siianoi  groups.  The  formation  of  amino  groups  may  also  cause 
accelerated  relaxation. 

Figure  3.5(b)  shows  that,  for  a  nitridation  temperature  of  800°C,  the  peak 
interfacial  nitrogen  concentration  does  not  appear  to  move  away  from  the 
interface  with  increasing  nitridation  time,  as  is  the  case  for  higher  mtridation 
temperatures.  Different  nitridation  kinetics  may  become  dominant  at  low 
mtridation  temperatures,  even  though  the  oxidation  temperature,  and  hence 
the  macroscopic  intrinsic  stress,  does  not  appear  to  affect  the  nitridation 
kinetics  which  determine  the  nitrogen  distribution.  The  fact  that  the  peak 
interfacial  nitrogen  concentration  is  observed  at  the  interface  itself  is  an 
indication  that  the  oxygen-rich  oxynitride  observed  at  the  interface  for  higher 
mtridation  temperatures  is  not  an  etching  artifact  (e.  g.  resulting  from 
preferential  etching  or  hydroxylation  of  the  surface  of  the  oxynitride  by  the 
etchant).  This  was  previously  confirmed  by  the  use  of  XeF2  gas  phase 
etching  of  Si  to  remove  the  substrate  and  expose  the  interface,13  and  by  the 
observation  of  the  relatively  oxygen-rich  interface  when  ion  beam  etching  is 
used.67'69 

III.  B.  Defect  Distributions 

Figure  3.6  shows  the  nitridation  time  (tN)  dependence  of  JAVFB|,  the  peak 
interfacial  nitrogen  concentration  [N]0,  and  the  distance  dNp  of  [N]p  from  the 
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Figure  3.6.  (a)  Magnitude  of  the  flatband  voltage  shift  |AVFJ  as  a  function  of 

nitridation  time  tN,  for  600-A-thick  oxides  nitrided  at  f000‘C  (from  Ref. 
57).  (b)  Peak  interfacial  nitrogen  concentration  [N]  and  distance  of 
[N]  from  the  interface  as  functions  of  tN,  for  a  100-A-tnick  oxide  nitrided 
at  iOOO'C  (from  Fig.  3.3). 
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interface.  Figure  3.6(a)  is  based  on  the  data57  of  Chen  et  at.,  taken  on 
600-A-thick  oxides  nitrided  at  1000‘C.  Unfortunately,  |A V FB |  was  not 
measured57  for  tN  <  30  minutes,  so  that  the  data  does  not  show  the  initial 
increase  in  |AVFB|,  as  is  observed  for  nitridation  temperatures  (TN)  belcw 
1000’C.  The  dashed  line  in  Fig.  3.6(a)  for  tN  <  30  minutes  is  the  expected 
behavior  of  |AVFBj,  based  on  the  initial  increase  in  (A VFB{  observed57  for  TN  < 
1000‘C.  The  time  scale  on  which  |AVFB|  increases  is,  however,  expected  to 
be  shorter  than  is  observed  for  TN  <  1000‘C  due  to  faster  reaction  rates  at 
elevated  temperatures.  Figure  3.6(b)  is  based  on  the  data  presented  in 
section  III. A  on  100-A-thick  oxides  nitrided  at  1000‘C.  Qualitatively  similar 
behavior  is  observed67"69  on  400-A-thick  oxides  nitrided  at  1100‘C.  Since 
the  data  represented  in  Figs.  3.6(a)  and  3.6(b)  were  taken  on  samples  with 
different  initial  oxide  thicknesses,  only  a  comparison  of  the  qualitative  trends 
is  useful.  In  the  initial  stages  of  nitridation,  [N]p  increases  rapidly  to  some 
limiting  value,  and  |AVfb|  increases.  As  tN  increases,  [N]  remains  nearly 
constant,  dNp  increases,  and  |AVFBj  increases  until,  after  some  value  of  dNo, 
i A V F B |  begins  to  decrease.  A  clear  correlation  between  the  nitrogen 
distnbution  and  |AVFB|  is  seen  to  exist. 

It  is  postulated  that  the  defects  giving  rise  to  the  charge  density 
responsible  for  | AV FB j  have  a  spatial  distribution  which  is  correlated  to, 
though  not  the  same  as,  the  nitrogen  distribution.  It  is  to  be  noted  that  the 
band  bending  represented  by  |AVFBj  depends  on  the  spatial  distribution  of  the 
charge  within  the  oxide.  The  oxide  charge  induces  an  image  charge  in  the 
substrate,  and  the  image  potential,  which  is  inversely  proportional  to  the 
charge  separation,  causes  the  band  bending.  Thus,  a  decrease  in  1 A VFB ] , 
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which  is  a  measure  of  the  band  bending,  could  be  caused  either  by  a 
decrease  in  the  total  charged  defect  density  or  by  an  increase  in  charge 
separation. 

The  prooosed  correlation  between  tne  charged  defect  cistriouticn 
responsible  for  i A V FB j  and  the  observed  nitrogen  distribution  may  be 
recognized  to  occur  as  follows.  As  the  substrate  reacts  with  the  diffusing 
species,  the  energy  of  the  system  is  lowered  by  the  formation  of  bonds 
between  Si  and  N  or  O.  The  presence  of  strain  at  the  SiC^/Si  interface, 
however,  partially  offsets  the  gain  in  energy  due  to  bond  formation.  This 
strain  energy  cost  can  be  lowered  if  defects  which  relieve  the  strain  locally 
are  formed  during  the  nitridation  process.  The  ability  to  lower  the  system 
energy  via  such  defects  thus  suggests  that  the  free  energy  of  defect  formation 
(AGf)  behaves  as  shown  schematically  in  the  inset  in  Fig.  3.7.  The  result 
would  be  an  increase  in  defect  concentration  at  the  interface  as  [N]p 
increases  in  the  initial  stages  of  nitridation.  As  the  nitridation  proceeds  and 
dNo  increases  by  virtue  of  the  movement  of  the  SiOxNy/Si  interface,  the 
defects  initially  formed  are  left  progressively  farther  from  the  interface.  A 
schematic  representation  of  this  model  is  shown  in  Fig.  3.7.  The  strain  energy 
released  via  defect  formation  is  then  expected  to  decrease  as  a  result  of  two 
different  effects:  (a)  the  decrease  in  strain  with  increasing  distance  from  the 
interface76'77  and  (b)  the  increasing  incorporation  of  nitrogen  in  the  film.  The 
significance  of  (b)  may  be  recognized  by  noting  that  chemical  vapor 
deposited  Sigf^  films  on  Si  are  known  to  be  in  macroscopic  intrinsic  tensile 
stress.82  Thus,  increasing  incorporation  of  nitrogen  in  the  S,02  film  may 
counter  the  compressive  stress  which  is  present  near  the  SiC^/Si  interface 
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Defect  Concentration  (arb.) 


Figure  3.7.  Schematic  representation  of  the  distribution  of  defects  at  various 
mtridation  times.  Inset:  Schematic  diagram  of  the  free  energy  of 
formation  AGf  vs.  configuration  coordinate  for  the  formation  of 
SiOxN  /Si  without  defects  (solid  line)  and  with  defects  (dashed  line)  to 
relievev  strain  energy. 
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prior  to  nitridation.  A  decrease  in  the  strain  implies  ihat  the  generation  of 
defects  is  no  longer  as  important  for  releasing  strain  energy,  ar,d 
consequently  the  defect  density  decreases,  as  shown  schematically  in  Fig. 

O  -* 
u.  /  . 

It  is  further  postulated  that  the  defects  discussed  above  act  as  charge 
traps.  This  postulate  has  the  following  consequences  for  the  electrical 
properties.  The  initial  increase  in  defect  density  results  in  an  increase  in  the 
charge  density,  a  concomitant  increase  in  band  bending,  and  hence  an 
increase  in  |AVFB|.  As  the  interface  moves  away  from  the  charged  defects,  the 
image  potential,  and  hence  |AVFB|,  decreases,  although  this  effect  could,  in 
principle,  be  offset  by  a  continuing  increase  in  the  charge  density.  As  the 
interface  continues  to  move  away  from  the  charged  defects,  however,  the 
decreased  defect  density  resulting  from  decreased  strain  will  cause  a 
decrease  in  the  charge  density  and  thus  a  decrease  in  |AVfb|.  In  addition  to 
the  effect  of  decreasing  strain,  electrical  properties  can  also  be  affected  by  a 
decreased  strain  gradient.  A  correlation  has  previously  been  noted83  84 
between  decreased  strain  gradients  and  lower  values  of  the  radiation-gener¬ 
ated  interface  trap  density.  A  decreased  strain  gradient  may  also  play  a  role 
in  the  behavior  of  |AVFB|. 

It  is  important  to  note  that  the  model  discussed  above  does  not  depend 
on  the  identity  of  the  defects  which  are  proposed  to  relieve  strain.  Some 
possibilities  are  dangling  orbitals  and  hydrogen-related  defects,  such  as  Si 
bonded  to  hydroxyl  or  amino  groups  or  to  hydrogen.  However,  if  the  defects 
are  hydrogen  related,  only  a  small  fraction  of  the  defects  could  be  electrically 
active,  since  the  concentration  of  hydrogen  in  nitrided  oxide  films  has  been 
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measured62  to  be  several  atomic  percent,  which  is  far  greater  than  the 
concentration  of  electrically  active  defects  required  to  account  for  the 
observed57  flatband  voltage  shifts. 

The  model  discussed  above  to  explain  the  behavior  of  IAVFQi  as  a 
function  of  nitridation  time  also  explains  the  behavior  of  |AVFB|  as  a  function  of 
other  processing  variables.  For  tN  =  30  minutes,  the  maximum  value  of  ;AVFB| 
occurs57  at  TN  =  1000‘C.  portN  =  60  minutes,  |AVFB|  reaches  a  maximum57  at 
Tn  =  950'C.  For  tN  =  90  minutes,  the  maximum  value  of  |AVFB|  still  occurs57  at 
Tn  =  950‘C  although  the  absolute  value  is  lower  than  for  tN  =  60  minutes.  For 
tN  =  1  SC  minutes,  the  maximum  value  of  |AVFB|  occurs57  at  TN  =  900‘C.  This 
behavior  is  a  consequence  of  slower  nitridation  reaction  rates  at  low 
temperatures.  This  can  be  s^un  by  noting  that  for  TN  =  800’C,  dNp  =  0  even 
for  tN  =  240  minutes,  as  shown  in  Fig.  3.5,  while  for  nitridation  temperatures  of 
1 000‘C  (see  Fig.  3.3),  1100’C, 67-69  and  1150’C  (see  Fig.  3.4),  dNp  becomes 
nonzero  within  tN  <15  minutes.  Thus,  at  lower  temperatures,  the 
nitridation-induced  defects  form  at  a  slower  rate  and  are  in  the  strained  region 
near  the  interface  longer,  resulting  in  |AVFB|  reaching  a  maximum  value  at  a 
later  time.  It  would  appear  that  a  nitridation  temperature  near  the  viscous  flow 
temperature  of  Si02  separates  fast  from  slow  nitridation  reaction  rates. 

It  has  also  been  reported85  that  under  the  same  nitridation  conditions, 
nitridation  of  thicker  oxides  results  in  increased  values  of  |AVFB|. 
Furthermore,  under  the  same  nitridation  conditions,  thicker  oxides  have  been 
reported64  to  have  less  nitrogen  incorporated  in  the  interfacial  region.  Results 
of  simulations  also  show  (see  section  IV. D)  a  reduced  interfacial  nitrogen 
concentration  for  a  larger  initial  oxide  thickness,  as  well  as  a  smaller  value  of 
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dND.  A  lower  value  of  [N]p  implies  a  lower  defect  concentration,  and  ti  ^  a 
lower  value  of  |uW FB j  within  the  modei  discussed  here.  However,  a  lower 
value  of  dNo  implies  a  higher  defect  concentration.  The  net  effect  may  be  an 
incrnased  defect  concentration  and  could  explain  the  observed35  larger 
values  of  iAVFBj. 

An  indication  of  the  defect  distribution  in  the  dielectric  film  can  also  be 
obtained  from  the  XPS  data  by  use  of  a  fluorine  marker.  The  HF/EtOH 
etchant  used  in  this  work  leaves  a  fluorine  residue  on  the  surface.  It  is 
suggested  that  the  fluorine  makes  bonds  preferentially  to  defects  in  the 
dielectric  film.  The  intensity  of  the  F  Is  photoelectron  peak  after  each  etch 
may  thus  yield  a  measure  of  the  defect  distribution  in  the  dielectric  film.  The 
surface  fluorine  concentration  should  be  considered  an  upper  limit  to  the 
defect  concentration,  since  some  of  the  fluorine  may  be  adsorbed,  possibly 
weakly,  at  sites  other  than  defect  sites. 

Figure  3.8  shows  the  number  of  F  atoms  bonded  to  the  dielectric  surface 
during  the  course  of  the  depth  profiling  of  samples  1  and  2,  corresponding  to 
the  data  in  Figs.  3.3(a)  and  3.3(b).  The  number  of  F  atoms  on  the  surface  was 
estimated  by  calculating  the  F  Is/Si  2 p  ratio  which  would  be  expected  for  a 
monolayer  (6.8  x  1014  atoms/cm2)  of  F  atoms  on  Si,  assuming  that  AF  (see 
Appendix  A)  is  1.  The  Si  2 p  intensity  is  estimated  by  calculating  the  expected 
attenuation  of  the  measured  clean  Si  substrate  signal  when  a  2.5  A  thick 
overlayer  with  a  25  A  escape  depth  is  present  on  the  surface.  The  absolute 
intensity  of  the  F  1  s  photoelectron  peak  for  a  monolayer  of  F  atoms  can  thus 
be  estimated  and  compared  to  the  actual  F  Is  intensities  measured  during 
the  course  of  the  depth  profile.  In  Fig.  3.8(a),  a  relatively  large  amount  of  F 
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Figure  3.8.  Surface  fluorine  concentrations  measured  during  the  course  of 
depth  profiles  of  samples  nitrided  at  1000'C  for  (a)  10  minutes,  and  (b) 
30  minutes. 


DISTANCE  FROM  SiOxNy/Si  INTERFACE  (A) 


bines  to  the  surface  when  the  interfacial  region  has  been  exposed  to  the 
etchant,  suggesting  the  formation  of  a  large  number  of  defects  at  the  interface 
in  the  initial  stages  of  nitridation.  Comparing  Fig.  3.8(b)  to  Fig.  3.8(b)  shows 
that  as  the  nitridation  proceeds  and  the  interfacial  nitrogen  peak  occurs  away 
from,  the  intedace,  the  maximum  surface  fluorine  concentration  is  lower  and 
also  occurs  further  from  the  interface.  This  is  precisely  the  behavior  expected 
of  the  defect  distribution  based  on  the  previously  proposed  strain-dependent 
defect  formation  model. 

Figure  3.9  shows  the  surface  fluorine  concentration  corresponding  to  the 
data  in  Fig.  3.4.  The  surface  fluorine  concentration  correlates  with  the 
nitrogen  distribution,  suggesting  an  increased  number  of  defects  as  the 
nitridation  time  increases  from  5  minutes  to  10  minutes,  but  distributed 
somewhat  further  from  the  interface.  An  increased  concentration  of  charged 
defects  would  cause  an  increase  in  the  tiatband  voltage  shift.  However,  an 
increase  in  the  separation  of  the  charge  from  the  interface  would  decrease 
the  flatband  voltage  shift.  Since  the  flatband  voltage  depends  on  both  the 
concentration  and  spatial  distribution  of  the  charged  defects,  it  is  not  clear 
how  the  electrical  measurements  would  be  affected.  The  electrical 
measurements57  on  which  the  defect  formation  model  is  based  do  not  cover 
the  nitridation  conditions  used  for  the  samples  corresponding  to  Figs.  3.4  and 
3.9. 

Figure  3.10  shows  the  surface  fluorine  concentration  corresponding  to 
the  data  in  Fig.  3.5.  Comparing  Figs.  3.10(a)  and  3.10(b)  to  Fig.  3.10(c) 
suggests  an  increased  defect  concentration  at  the  interface  as  the  nitridation 
time  increases  from  1  hour  to  4  hours.  This  is  consistent  with  the  reported5' 
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Figure  3.9.  Surface  fluorine  concentrations  measured  during  the  course  of 
depth  profiles  of  samples  nitrided  at  1 1 50 ‘C  for  (a)  5  minutes,  and  (b) 
10  minutes. 
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Figure  3.10.  Surface  fluorine  concentrations  measured  during  the  course  of 
deoth  profiles  of  samples  (a)  oxidized  at  825*C,  then  mtrided  at  800'C 
for  60  minutes,  (b)  oxidized  at  1000’C,  then  mtrided  at  800'C  for  60 
minutes,  and  (c)  oxidized  at  1000’C,  then  nitrided  at  800’C  for  240 
minutes. 


r.ionotonic  increase  in  the  magnitude  of  the  fiatband  voltage  shift  with 
increasing  mtndation  time  at  low  nitridation  temperatures.  Comparing  Figs. 
3.10(a)  anc  3.10(b)  to  Fig.  3.5(a)  suggests  that,  even  though  the  oxidation 
temperature  aoes  not  appear  to  affect  the  nitridation  Kinetics  responsible  for 
the  nitrogen  distribution,  it  does  affect  the  defect  formation  kinetics.  The  data 
in  Fig.  3.10  suggest  that  defect  formation  is  reduced  in  the  oxide  which  is 
expected  to  have  the  higher  strain.  It  has  previously  been  suggested83'84 
that  increased  strain  in  the  bulk  of  the  oxide  reduces  the  interfacial  strain 
gradient,  resulting  in  a  reduced  concentration  of  radiation-generated  interface 
traps.  An  analogous  effect  may  be  occurring  in  the  present  work. 

Figures  3.8  through  3.10  show  that  the  surface  fluorine  concentration  in 
mtrided  oxides  is  near  1014crrr2.  This  is  far  greater  than  the  number  of 

charged  defects  required  to  account  for  observed57  fiatband  voltage  shifts. 

* 

For  comparison,  the  surface  fluorine  concentration  observed  during  the 
course  of  a  depth  profile  of  a  thermally  grown  SiC>2  film  on  Si  is  near  1013 
cm'2,  and  is  near  1012  cm'2  on  the  surface  of  the  clean  Si  substrate  after  the 
S1O2  has  been  completely  removed.  However,  device  quality  thermal  oxides 
have  a  concentration  of  electrically  active  defects  as  low  as  1010  eV'1  cm'2  at 
midgap.  It  is  clear  that  only  a  small  fraction  of  the  defects  associated  with  the 
fluorine  marker  can  be  electrically  active.  It  is  suggested,  however,  that  the 
distribution  of  electrically  active  defects  may  follow  the  same  behavior  as  the 
distribution  of  defects  deduced  from  the  surface  fluorine  concentration. 
Comparing  Figs.  3.8  to  3.10  with  Figs.  3.3  to  3.5,  it  is  clear  that  the  surface 
fluorine  concentration,  and  hence  the  defect  distribution,  is  correlated  to  the 
nitrogen  distribution.  This  is  consistent  with  the  strain-dependent  energy  of 


formation  of  defects,  proposed  earlier  to  explain  the  observed57  behavior  of 
the  flatband  voltage  shift.  In  addition,  recent  electrical  measurements58  show 
that  the  positive  charge  is  centered  24  A  from  the  SiOxNy/Si  interface  for  a 
100  A  Si02  film  which  has  been  nitnded  at  1000°C  for  10  minutes.  The  fact 
that  the  peak  interfacial  nitrogen  concentration,  the  maximum  surface  fluorine 
concentration,  and  the  positive  charge  center  all  occur  away  from  the 
interface  strongly  suggests  that  the  defects  are  associated  with  the  nitrogen 
and  that  the  fluorine  makes  bonds  preferentially  with  defect  sites,  thus 

i 

supporting  the  previously  proposed  defect  model. 

IV.  Modeling  of  the  Kinetics  of  Nitridation 

This  section  describes  simulations  of  the  nitridation  process.  The 
nitridation  process  is  described  in  section  IV. A,  the  kinetic  model  presented  in 
section  IV. B,  the  choice  of  kinetic  parameters  described  in  section  IV. C,  anc 
the  simulation  results  are  presented  in  section  IV. D. 

IV.  A.  The  Nitridation  Process 

Several  reactions  are  involved  in  the  nitridation  process.  These  include 
reaction  of  the  Si02  with  NH3,  diffusion  of  the  NH3  and  reaction  byproducts 
through  the  oxynitride,  possible  reaction  of  the  byproducts  with  the  SiOxNy, 
and  reaction  of  the  substrate  with  the  diffusing  species.  In  the  following,  the 
nature  of  each  of  these  reactions  will  be  considered. 

There  are  three  stable  phases  in  the  Si  -  N  -  O  phase  diagram86  --  S1O2. 
S12N2O,  and  Si3N4.  The  reaction  of  S1O2  with  NH3  may  therefore  be 
described  by  the  following  reactions: 


2  Si02  +  2  NH3  =  Si2N20  +  3  H20, 

3  Si02  +  4  NH3  =  Si3N4  +  6  H20. 

A  recent  infrared  (IR)  study61  found  that  nitrided  oxides  have  an  IR  absorption 
at  a  frequency  which  is  near  that  observed  for  trisilyiamine  (N(SiH3)3).  The 
bonding  about  the  nitrogen  in  trisilyiamine  is  trigonal  planar,  rather  than  the 

more  common  pyramidal.  The  IR  study61  would  not  distinguish  between 

! 

Si2N20  and  Si3N4,  however,  since  the  Si-N-Si  bond  angle  in  Si2N20  is3'88 
120*.  while  that  in  Si3N4  is88  119.9',  both  consistent  with  trigonal  planar 
bonding.  In  the  depth  profiles  of  nitrided  oxides  presented  in  Section  III. A,  the 
concentration  ratio  [N]/([0]  +  [N])  never  exceeds  0.67,  consistent  with  the  final 
reaction  product  being  Si2N20.  Of  the  two  reactions  being  considered,  (3.1) 
is  thermodynamically  favored.60'65  Finally,  Si2N20  has  been  reported89  to 
form  when  Si02  is  heated  in  NH3  in  the  temperature  range  550-1 250’C, 
while  formation  of  Si3N4  is  reported89  at  higher  temperatures.  We  therefore 
conclude  that  the  formation  of  Si2N20  via  reaction  (3.1)  is  more  likely  to 
occur  for  the  nitridation  conditions  being  considered  here. 

The  diffusion  of  gases  in  fused  silica  has  been  extensively  studied.71'73 
H20  reacts  with  Si02  to  form  silanol  via  the  reaction 

Si  -  O  -  Si  +  H20  =  2  Si  -  OH.  (3.3) 

This  reaction  is  reversible,  with  the  H20  being  the  diffusing  species.  A  similar 
reaction  presumably  occurs  with  NH3: 


(3.1) 

(3.2) 
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This  differs  from  the  case  of  H2O,  however,  since  the  amino  groups  may 
become  immobile  by  the  formation  of  multiple  Si  -  N  bonds,  without 
necessarily  releasing  the  diffusing  NH3>  via  the  reactions 

Si  -  NH2  +  Si  -  0  -  Si  =  Si2  -  NH  +  Si  -  OH.  ( 3 . 5 1 

Si-NH2  +  Si-0H  =  Si2-NH  +  H20,  (3.6) 

2  Si  -  NH2  =  Si2  -  NH  +  NH3  (3.7) 

and  similar  reactions  involving  the  formation  of  a  third  Si  -  N  bond.  These 

reactions  may  be  reversible.  However,  the  necessity  of  breaking  multiple  Si  - 
N  bonds  before  diffusion  can  occur  would  be  expected  to  greatly  reduce  the 
rate  of  formation  of  the  diffusing  NH3,  compared  to  the  rate  in  the  case  of  Si  - 
NH2,  which  requires  breaking  a  single  Si  -  N  bond. 

The  occurrence  of  reactions  such  as  (3.6)  means  that  H20,  as  well  as 
Nh  o,  is  present  and  diffusing  in  the  SiOxNy  during  the  nitridation.  The 
presence  of  both  diffusing  species  makes  it  unlikely  that  Si02  or  Si3N4  will 
result  from  the  reaction  of  the  substrate  with  the  diffusing  species.  The 
formation  of  Si2N20  could  result  from  reaction  of  the  substrate  directly  with 
the  diffusing  species,  via  the  reaction60’65 


2  Si  +  2  NH3  +  H20  =  Si2N20  +  4  H2. 


(3.8) 


Alternatively,  the  interfacial  Si02  may  also  be  involved,  via  the  reaction65 


3  Si  +  Si02 +4  NH3  =  2  Si2N20 +  6  H2.  (3.9) 

Of  these  two  reactions,  (3.8)  is  energetically  more  favorable.65  Depending  on 
the  relative  concentrations  of  the  diffusing  species,  nowever,  oxynitrides  of 
different  average  stoichiometries  could  result. 

The  nitridation  process  can  be  summarized  as  follows.  The  ammonia 
enters  the  Si02,  through  which  it  can  diffuse  and  chemically  react.  The 
substitution  of  N  for  0  in  the  oxide  can  release  H20,  which  is  then  free  to 
dinuse  and  possibly  react  with  the  oxynitride.  Those  diffusing  species  which 
do  not  become  immobilized  by  reaction  with  the  oxynitride  and  which  diffuse 
to  the  interface  can  react  with  the  substrate.  A  kinetic  model  of  the  various 
reactions  discussed  above  which  are  involved  in  the  nitridation  process  will 
now  be  described. 

IV.  B.  The  Kinetic  Model 

The  processes  discussed  in  Section  IV. A  can  be  modeled  using  either  a 
numerical  approach  or  by  attempting  an  analytical  solution  of  the  diffusion 
equation  with  chemical  reaction.  The  diffusion  equation  with  chemical 
reaction  has  been  solved  for  several  cases.90  However,  the  cases 
considered90  involved  a  constant  diffusion  coefficient.  In  the  present  case, 
the  diffusion  coefficient  will  depend  upon  the  amount  of  nitrogen  incorporated 
in  the  oxynitride,  since  nitrided  oxides  are  known50'52-53'60  to  be  better 
diffusion  barriers  than  is  Si02.  From  the  data  presented  in  Section  III. A.  it  is 


48 


clear  that  the  nitrogen  concentration  varies  both  with  the  nitridation  time  and 
the  position  within  the  oxide.  The  diffusion  coefficient  will  therefore  also  be  a 
function  of  position  and  time,  greatiy  complicating  the  solution  even  for 
simoiified  cases.  Adding  the  effects  of  strain  and  multiple  diffusing  species 
(e.g.  to  account  for  reaction  (3.6))  further  complicates  the  equation  to  De 
solved.  A  numerical  method  then  becomes  preferable  for  reasons  of 
simplicity  and  flexibility. 

Since  amorphous  materials  (SiC^,  SiOxNy)  are  being  considered,  a 
model  which  considers  local  atomic  arrangements  is  impractical,  especially 
for  the  size  ('•100  A  thick  SiC>2  layer)  of  the  system  being  modeled.  It  is 
therefore  necessary  to  consider  only  average  properties.  This  makes  dealing 
with  a  large  number  of  each  reacting  species  feasible,  which  is  necessary  if 
the  probabilities  fui  the  various  reactions  (diffusion,  incorporation)  vary 
widely.  In  addition,  a  model  which  considers  only  average  properties  and  a 
large  number  of  each  reacting  species  will  result  in  less  statistical  scatter  in 
the  final  results  than  a  Monte-Carlo  approach  in  which  every  atom  is 
considered. 

Assume  that  the  SiOxNy  is  laterally  uniform  both  chemically  and 
structurally  at  all  times.  The  dielectric  layer  can  then  be  divided  into  a  number 
of  thinner  layers,  each  with  its  own  average  diffusion  coefficients, 
incorporation  probabilities,  and  chemical  composition.  This  reduces  the 
problem  to  one  dimension.  We  shall  assume  that  the  hopping  distance  for 
diffusing  species  is  2.63  A,  which  is  the  O  -  0  nearest  neighbor  separation  in 
a-quartz.  Letting  the  layer  thickness  be  the  same  as  the  hopping  distance,  35 
layers  would  then  correspond  to  a  92-A-thick  dielectric  layer. 
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Before  discussing  the  simulated  nitridation  in  detail,  an  overview  is 
presented  in  the  following.  After  each  time  interval  At.  the  amount  of  diffusing 
nitrogen  in  the  surface  layer  is  initialized  to  the  solubility  limit.  Each  diffusing 
SDecies  (NH3  and  H2O)  is  then  allowed  a  single  incorporation  attempt  in  the 
oxynitnae  or  substrate  layer  in  whicn  it  is  located.  That  fraction  of  each 
diffusing  species  which  does  not  react  is  then  allowed  a  single  interlayer 
hopping  attempt,  and  the  process  repeats.  An  overview  of  the  kinetic 
simulation  is  shown  schematically  in  Fig.  3.11.  The  probability  of  each  of 
these  reactions  occurring  within  At  can  be  calculated,  as  discussed  later.  In 
the  following,  the  probabilities  will  be  assumed  to  be  known. 

Let  N0  be  the  total  number  of  available  hopping  sites  within  each  layer, 
and  f0  the  fraction  of  available  sites  in  the  surface  layer  which  are  occupied  by 
the  diffusing  N  species  at  time  t  =  t0,  when  nitridation  is  initiated.  Let  Nh  (t) 
and  Ninc  (t)  be  the  amount  of  nitrogen  hopping  and  incorporated,  respectively, 
within  the  j,h  layer  at  time  t,  where  j  =  1  is  the  surface  layer.  Nhop  1  (t0)  is  then 
equal  to  f0N0  and  Njnc  (t)  is  zero  for  all  j.  Since  the  surface  layer  is  in  direct 
contact  with  the  NH3  ambient,  f0  will  be  chosen  so  that  f0N0  is  the  solubility 
limit  of  the  diffusing  N  species  in  SiC^,  as  discussed  in  the  next  section.  As 
the  nitridation  proceeds,  the  incorporation  of  nitrogen  will  decrease  the 
solubility  limit,  since  nitrided  oxides  are  known  to  be  good  diffusion 
barriers. 50'52'53'60  It  is  assumed  that  the  solubility  limit  for  the  diffusing  N 
species  decreases  linearly  as  Njncl(t)  increases  to  some  limiting  value,  i.e. 
the  final  reaction  product  is  assumed  to  be  a  perfect  diffusion  barrier.  The 
effect  of  this  assumption  is  to  halt  the  nitridation  as  Njnc1(t)  asymptotically 
approaches  the  limiting  value.  The  choice  of  this  limiting  value  determines 
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the  final  reaction  product  in  tne  simulated  nitridation  process.  The  limiting 
value  will  be  written  as  xLN0,  so  that  the  final  reaction  product,  Si3N4x06(1.x),  is 
obtained  when  x  =  xL.  The  amount  of  diffusing  nitrogen  in  the  surface  layer  is 
then  determined  by 

NhoB.i  (')  -  fo<No  -  Nmc.i  (3-10: 

The  process  of  incorporation  is  shown  schematically  in  the  portion  of  Fig. 
3.11  enclosed  by  the  dashed  lines.  The  amount  of  diffusing  nitrogen  which 
incorporates  in  the  jth  oxynitride  layer  within  the  nth  time  interval  At  is  given  by 

ANinc,i(tn)  =  Nhop)(tn)PNoxjnCiJ(At),  (3.1  r 

where  tn  =  t0  +  nAt  and  PN  QX  incj(At)  is  the  probability  of  nitrogen  incorporation 
m  the  j,h  oxynitride  layer  within  time  At.  The  total  amount  of  nitrogen 
incorporated  in  the  j,h  oxynitride  layer  is  then 

^inc.j^n^  =  ^Njncj(tn)  +  NinCtj(tn_.,).  (3.12) 

Nhopj(tn)  is  correspondingly  decreased  by  ANjncj(tn).  From  reactions  (3.1 )  and 
(3.2),  each  nitrogen  which  incorporates  in  the  dielectric  layer  displaces  on  the 
average  3/2  O.  Therefore,  ^hop,  ■  (tn)  is  increased  by  3/2  ANjnC](tn). 
Incorporation  of  oxygen  in  the  oxynitride  is  similar. 

Incorporation  of  nitrogen  and  oxygen  in  the  substrate  is  only  slightly 
more  involved  (see  Fkw  3.11).  For  a  substrate  layer  which  has  only  partially 
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been  converted  to  an  oxynitride,  the  probability  that  the  diffusing  species 
occupy  sites  which  have  already  been  oxidized  or  nitrided  is  calculated.  That 
fraction  of  each  diffusing  species  contained  in  oxynitride  areas  is  treated  in 
exactly  the  same  manner  as  discussed  above.  The  fraction  which  occupies 
unreacted  substrate  areas  is  treated  in  a  similar  manner,  using  appropriately 
different  probabilities  and  with  no  displacement  to  create  diffusing  species. 

The  fraction  of  each  diffusing  species  which  has  not  incorporated  is  next 
allowed  to  diffuse,  as  shown  schematically  in  the  portion  of  Fig.  3.1 1  enclosed 
Dy  oot-dashed  I  nes.  In  this  model,  only  single-layer  hops  are  allowed,  and 
only  one  interlayer  hop  is  allowed  within  the  time  interval  At.  The  probability 
that  a  hop  will  occur  is  decreased  by  the  probability  of  target  sites  alreacy 
being  occupied  by  diffusing  species,  and  by  the  fact  that  some  of  the  target 
sites  are  part  of  a  diffusion  barrier.  This  model  results  in  a  diffusion  constant 
which  decreases  linearly  with  increasing  N  incorporation  until  it  is  zero  for  x  = 
xL.  It  is  assumed  that  the  diffusing  species  and  the  SiOxNy  form  an  ideal 
solution,  with  no  heat  of  mixing  and  no  change  in  chemical  potential  as  a 
function  of  position  in  a  concentration  gradient,  so  that  a  chemical  potential 
does  not  provide  a  driving  force  for  diffusion.  The  simplifying  assumptions  of 
no  interaction  between  diffusing  species  and  no  dependence  of  the  diffusion 
constant  on  the  concentration  of  diffusing  species  are  also  made.  The 
amount  of  nitrogen  which  hops  from  layer  j  in  the  +  (towards  the  substrate) 
and  -  (towards  the  surface)  directions  is  then  given  by 

AIW,(t„)  =  (1/2)Nho01(yPNhopj51(tn){l  - 
(N,„c,i=I<g'*L  -  Nh0D ij±1  (t>  Ohop  (tJJ/Nj.  (3. 1 3) 
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The  factor  of  1/2  gives  equal  weight  to  forward  and  backward  hopping,  i.  e. 
only  self-diffusion  is  considered,  and  there  are  no  driving  forces  which  result 
in  a  Dreferent.al  hopDing  direction.  Nh0D  ;(tn)  is  decreased  and  Nh0D  |r,(tri)  are 
increased  Dy  ANhop  ),i(tri).  Similar  equations  determine  oxygen  hopping. 
Diffusing  species  which  occupy  unreacted  substrate  sites  are  not  allowed  to 
hop  further  into  the  substrate,  i.e.  the  solubility  of  the  diffusing  species  in  the 
substrate  is  assumed  to  be  negligible. 

Each  of  the  probabilities  mentioned  above  must  be  estimated.  In  order  to 
calculate  the  probability,  it  is  assumed  that  each  of  the  reactions  considered 
follows  first  order  chemical  kinetics.  The  concentration  of  a  reacting  species 
as  a  function  of  time  c(tn)  is  then  given  by 

c(tn)  =c(tn.1)exp(-RAt),  (3.14) 

wnere  R  is  the  rate  constant.  W'thin  time  At,  the  probability  P(At)  for  the 
reaction  to  occur  is 


P(At)  =  1  -  c(\nmnA)  =  1  -  exp(-RAt).  (3.15) 

Assuming  an  Arrhenius  form  for  the  temperature  dependence  of  the  reaction 
rate,  Eq.  (3.15)  becomes 

P(At)  =  1  -  exp[-AtR0exp(-AEA/kBT)j,  (3.16' 
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where  aEa  is  the  activation  energy  and  kB  is  Boltzman's  constant.  Thus,  if 
AEa  and  R0  can  be  determined,  the  probability  of  a  reaction  occurring  within 
time  At  can  be  calculated.  The  choice  of  these  kinetic  parameters  will  now  be 
consiaered. 

IV.  C.  Choice  of  Kinetic  Parameters 

The  solubility  limits  in  fused  silica  of  gases  which  do  not  readily  react 
with  the  SiOg  matrix,  such  as  He,  Og.  and  N2,  are  generally  0.01  -  0.03  times 
the  concentration  of  the  atoms  or  molecules  in  the  gas  phase.'1  For  reactive 
H2O.  however,  the  solubility  limit  is  three  orders  of  magnitude  higher71'73  due 
to  reaction  (3.3).  It  is  estimated72  that  the  solubility  of  free  H2O  is  about  the 
same  as  other  gases  in  Si02<  with  the  apparent  higher  solubility  due  to  the 
fact  that  most  of  the  H2O  is  trapped  as  silanol  at  any  given  time.  Experiments 
have  shown73  that  the  solubility  of  H2O  in  fused  silica  is  nearly  constant  in  the 
temperature  range  900-1 200'C,  and  increases  at  lower  temperatures. 
Expressed  as  the  number  of  hydroxyl  groups  per  Si  atom  in  SiC>2,  the 
solubility  limit  of  H2O  ranges73  from  3x1  O'3  at  1200'C  to  4x1 0'3  at  800'C, 
which  is  the  mtridation  temperature  range  of  interest  in  this  work.  The 
ocr,,rrence  of  reaction  (3.4)  is  likely  to  make  the  solubility  of  NH3  in  SiC>2 
close  to  that  of  H2O,  or  possibly  somewhat  lower  due  to  the  larger  size  of 
NH3  compared  to  H2O.  For  these  simulations,  f0  was  chosen  to  be  2x1  O'3  at 
high  mtridation  temperatures,  and  2.7x1  O'3  below  900‘C.  Since  f0  is  defined 
relative  to  the  number  of  nitrogen  incorporation  sites,  and  4  nitrogen  atoms 
can  incorporate  for  every  3  Si  atoms,  the  corresponding  number  of  amino 
groups  per  Si  atom  is  2.7x1  O'3  for  mtridation  temperatures  of  900'C  and 
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above,  and  3.6x1  O'3  at  lower  temperatures.  For  reasons  discussed  in  Section 
IV. A,  the  finai  nitridation  reaction  product  is  likely  to  be  S^^O,  so  xL  is  set  to 
0.75  in  the  simulation.  The  value  of  At  is  chosen  to  be  large  enough  to 
minimize  computing  time,  yet  small  enough  to  prevent  the  reaction 
prooaoiiities  from  approaching  1.  The  value  of  At  is  therefore  chosen  to  be  in 
the  range  i0'2  -  10'4  seconds,  depending  on  the  nitridation  temperature.  Tne 
value  of  N0  is  chosen  to  be  1012,  which  would  correspond  to  a  surface  area  of 
the  order  of  1 0'4  cm2. 

To  the  author's  knowledge,  no  experimental  measurements  of  the  kinetic 
parameters  for  reactions  (3.1)  and  (3.2)  have  been  reported.  R0  and  AEA  for 
nitrogen  incorporation  in  the  oxide,  R0_(N ,inc.ox)  and  AEA  (N  inc _0X)I  can  be 
estimated  from  measured  surface  nitrogen  concentrations  from  samples 
mtrided  at  1000'C  (Fig.  3.3)  and  800‘C  (Fig.  3.4)  for  times  sufficier'y  short 
that  the  surface  has  not  reached  its  saturation  nitrogen  concentration. 
Equation  (3.16)  can  be  written  as 


R0(T1)  =  -(1/At)exp[AEA/kaT.]ln[1  -  P(At,  T^].  (3.17 

A  similar  equation  can  be  written  at  T  =  T2.  Assuming  that  R0  is  independent 
of  temperature,  the  two  equations  can  be  set  equal  to  each  other  to  obtain 


AEa  =  [KgTJ^  -  T1  )]ln[P(At,  T2)/P(At.  T1 )],  (3. 1 S 


which  is  valid  if  the  probabilities  are  much  less  than  1  Estimates  for  the 
probabilities  at  T.  and  T2  are  obtained  by  using  the  assumed  values  of  f0  and 
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adjusting  the  probabilities  to  obtain  the  observed  surface  nitrogen 
concentrations  at  the  temperatures  and  times  considered.  This  gives 

^A,N,nc.ox)  =  '-47  eV  and'  from  EP-  (3-1 7),  R0.(N,nc.ox)  =  3.65x1 05  sec'1.  If  the 
temperature  dependence  of  the  solubility  is  taken  into  account,  AEa  fN  nc  ox;  = 

1 .6  eV  ano  R0  (N  incox)  =  1 ,2x106  sec-'  are  obtained. 

Estimates  of  the  kinetic  parameters  for  the  reverse  reaction  are  more 
difficult  to  obtain.  For  each  Si-0  bond  which  is  replaced  by  a  Si-N  bona  via 
reaction  (3.1)  or  (3.2).  the  free  energy  is  increased  by  -C.2  eVA0-53  If  the 
energy  barrier  is  unchanged  by  the  nitridation  reaction.  AEa  for  the  reverse 
reaction.  AEAi0incox),  would  therefore  be  0.2  eV  less  than  AEA  ;N,ncox  . 
Relaxation  of  the  local  atomic  environment  to  a  more  stable  arrangement  after 
the  nitridation,  however,  may  affect  the  barrier  height.  For  those  calculations 
considering  a  reversible  nitridation,  R0(Oincox)  is  assumed  to  be  the  same  as 

^O.iN.inc.ox)*  and  dlEA.(0.mc.ox)  'S  ^  '4  e'7- 

Macroscopic  stress  is  assumed  not  to  affect  the  incorporation  reactions 
discussed  above,  but  the  effects  of  interfacial  strain  are  included  in  the 
simulation.  Since  the  strained  bonds  are  more  easily  broken, 7S’77-:’9-80  a 
lower  activation  energy  for  nitrogen  incorporation  in  the  interfacial  region  is 
assumed,  decreasing  linearly  by  0.3  eV  as  the  interface  is  approached  in  the 
10  dielectric  layers  (26.3  A)  closest  to  the  interface,  to  simulate  the  bond 
strain  gradient.75'77  Enhanced  nitrogen  incorporation  in  the  i nterfacial  region 
may  be  an  indirect  consequence  of  strain,  since  increased  hyarogen 
incorooration  near  the  interface64-91  may  act  as  a  strain  relief  mechanism. 
Reclacmg  Si-OH  or  Si- H  bonds  with  Si-N  bonds  may  be  energetically  more 
‘avorable  than  reaction  (3.4). 


The  activation  energy  for  the  reaction  of  the  diffusing  H2O  with  the 
suDstrate,  AE,  ..n  can  be  obtained  from  kinetic  studies  of  the  steam 
oxidation  of  Si  (1 00).  The  oxide  thickness  increases  linearly92  with  time  in  the 
thin  oxide  regime  where  the  oxidation  rate  is  limited  by  the  interfacial 
reaction,  and  parabolically92  for  thicker  oxides  where  the  oxidation  is  diffusion 
limited.  In  the  regime  where  the  oxidation  is  reaction  limited,  the  activation 
energy  has  been  measured93  to  be  2.05  eV  for  oxidation  temperatures  above 
900‘C  and  1.6  eV  for  oxidation  temperatures  of  900'C  and  below.  A  similar 
temperature  effect  has  been  observed94'95  for  oxidation  of  Si  in  O2  and  has 
been  attributed  to  the  macroscopic  stress  observed74  in  low  temperature 
oxides.  In  the  simulation,  R0  (o  me  si)  is  set  arbitrary  lar9e  PO12  sec'1)  to 
ensure  that  the  reaction  is  diffusion  limited.  Although  an  oxide  thickness  of 
100  A  is  not  in  the  diffusion  limited  regime,  this  choice  is  motivated  by  the  fact 
that  nitrided  oxides  are  good  diffusion  barriers,50'52'53'60  and  nitrogen 
incorporates  rapidly  in  the  interfacial  region,13'1 8'67'69  thus  forming  a  diffusion 
barrier  early  in  the  nitridation  process. 

The  kinetics  of  the  reaction  between  Si  and  NH3  are  poorly  understood, 
though  some  studies  have  been  done.59'60'64'67'96'97  Attempts  to  determine 
the  linear  rate  constant,  which  is  determined  by  the  interfacial  reaction,  are 
complicated  by  the  fact  that  the  thickness  of  the  grown  Sig^  film  is 
self-limiting  at  a  tew  tens  of  Angstroms,  and  the  grown  films  have  a  significant 
oxygen  content  which  may  affect  the  nitridation  kinetics.  In  the  simulations, 
the  kinetic  parameters  for  the  reaction  of  the  substrate  with  NH3,  R0  ,N  mc  Sii 
and  AEi  ,N  ,nc  Si|,  are  chosen  to  be  the  same  as  for  the  reaction  with  H20. 

As  oreviously  mentioned,  the  diffusion  of  HoO  in  fused  silica  has  been 
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extensively  studied.71'73  Reaction  (3.3)  is  not  a  first  order  reaction,  but  rather 
depends  on  the  square  root  of  the  H2O  concentration,  as  expected  if  reaction 
(3.3)  is  an  elementary  reaction.  Reaction  (3.3)  is  not  the  only  pathway  by 
which  H9O  diffusion  can  occur,  however.  Reactions  such  as  (3.6)  can  aisc 
occur,  and  hydroxyl  groups  in  a  strain  gradient  can  diffuse  singly/'  rather 
than  in  pairs  as  expected  from  reaction  (3.3).  The  existence  of  multiple 
reaction  pathways  makes  estimating  the  overall  order  of  the  reaction  difficult. 
Fcr  simplicity,  first  order  kinetics  will  be  assumed. 

The  activation  energy  for  diffusion  of  H2O  in  fused  silica  is  0.78  eV,  in 
exact  agreement  with  measurements93  of  the  parabolic  (diffusion  limited)  rate 
constant  for  oxidation  of  Si  in  H2O  above  900‘C.  At  900'C  and  below,  the 
activation  energy  is93  1.17  eV.  In  the  simulations,  the  activation  energy  for 
diffusion  of  H2O  in  the  oxynitride,  AEa  ;0hopox),  is  set  t0  these  values  in  the 
appropriate  temperature  ranges.  R0.(O.hop.ox)  can  be  estimated  from  the 
effective  diffusion  constant  Deff,  which  is157  6.3x1  O'10  cm2  sec'1  at  1000'C. 
using  the  relation98 


Deff  =  (1/2)Ra2,  (3.19) 

where  a  is  the  hopping  distance.  Using  the  Arrhenius  relationship  between  R 
and  R0  and  AEA,,0hODOX;  =  0.78  eV,  which  is  appropriate  at  1000'C. 
R 3  ' o  hco  ox  =  2.2x1 09  sec'1  is  obtained.  In  the  simulations.  Ro.^o.noo.ox)  IS 
assumed  to  be  independent  of  temperature,  and  only  AEa  (Q  hopox)  is  changed 
at  low  mtridation  temperatures. 

Interfacial  strain  may  also  affect  the  diffusion  kinetics.  The  smaller  ring 


size  m  the  strained  region  should  increase  the  activation  energy  for  diffusion 
of  species  which  do  not  react  with  the  SiC>2  matrix.  In  fact,  the  kinetics  of  the 
dry  O2  oxidation  of  Si  has  been  explained99'100  with  an  interfacial  layer  in 
which  a  lower  diffusivity  is  exhibited,  and  a  higher  activation  energy  uas  been 
measured'01  for  thin  oxides,  Enhanced  diffusivity  has  been  proposed,7- 
however,  for  hydroxyl  groups  in  the  bond  strain  gradient  model.  Enhanced 
diffusivity  m  thin  (<  100  A)  oxides  via  transport  through  micropores  has  also 
been  proposed102  to  explain  the  initial  oxidation  of  Si  in  O2.  which  is  faster 

1 

than  the  linear-paraboiic  model  predicts.  Increased  incorporation  of 
hydrogen  near  the  interface64-91  may  contribute  to  a  more  open  oxiae 
structure  in  which  diffusion  is  enhanced.  The  effect  of  mterfaciai  strain  on 
diffusion  is  therefore  simulated  by  linearly  decreasing  the  activation  energy 
(i.e.  increasing  the  hopping  probability)  by  0.15  eV  as  the  interface  is 
approached  in  the  10  dielectric  layers  closest  to  the  interface. 

The  activation  energy  for  ammonia  diffusion  in  fused  silica,  AEa  (N  hop  ox), 
has  not  been  reDorted,  but  can  be  estimated.  Activation  energies  for  diffusion 
in  fused  silica  for  many  gases  follow  the  relation71 

(AEA)1/2  =  a  +  br,  (3.20' 

wnere  a  and  b  are  constants  and  r  is  the  molecular  radius.  However, 
molecular  size  alone  does  not  determine  the  activation  energy.  For  example. 
H20  (r  =  1 .65  A,  AEA  =  0.78  eV)  is  close  in  size71-72  to  Ar  (r  =  1 .6  A,  AEa  = 
1.16  eV),  02  (r  =  1.6  A,  AEa  =  1.09  eV),  and  N2  (r  =  1.7  A,  AEa  =  1 .12  eV),  but 
the  activation  energy  is  lower  than  expected  from  size  considerations  alone. 


Pari  of  the  problem  may  be  the  measurement  of  the  molecular  radius,  which 
is  typically  calculated  from  gas  viscosities  assuming  a  spherical  molecule, 
which  is  inappropriate  for  HpO. 

An  alternative  explanation  is  that,  of  the  gases  studied,71-72  only  l-UO  nas 
a  dipole  moment  (1.85  debyej.  The  interaction  of  the  dipole  with  the 
polarizable  S1O2  matrix  could  lower  the  activation  energy.  The  activation 
energy  for  NH3  diffusion  in  SiC>2  based  on  molecular  size  (r  =  1.82  A)  alone 
would  be  ~1.5  eV.  However,  since  NH3  has  a  dipole  moment  (1.47  debye) 
and  is  likely  to  diffuse  by  a  mechanism  similar  to  that  of  H2O,  a  lower 
activation  energy  would  be  expected. 

In  the  simulations,  AEA;Nhopox)  is  set  to  1.2  eV  for  nitridation 
temperatures  above  900'C,  while  1.4  eV  is  used  for  nitridation  temperatures 
of  900’C  and  below.  R0(Nhopox)  is  adjusted  to  reasonably  reproduce  data 
from  1 000’C  nitridations,  then  fixed  at  the  same  value  at  other  temperatures. 
In  tne  simulations,  R0{Nhopox)  =  3. 5x1 07  sec1  has  been  used.  This  gives  a 
diffusion  coefficient  of  -2x1  O'13  cm2  sec'1  at  1000’C,  considerably  larger  than 
the  previously  reported65  10'16-10'15  cm2  sec'1.  However,  the  mathematical 
treatment  used  in  that  study65-90  neglects  the  chemical  reaction  occurring  and 
treats  the  measured  concentration  of  nitrogen  incorporated  in  the  oxynitride 
as  the  concentration  of  the  diffusing  species.  For  reasons  previously 
discussed,  however,  the  concentration  of  diftusing  NH3  is  likely  to  be  quite 
small  (f0  =  2x1 0‘3),  and  the  measured  incorporated  nitrogen  concentration 
consists  primarily  of  the  immobile  nitrogen  which  has  formed  multiple  Si-N 
bonds.  In  the  simulations,  the  effect  of  interfacial  strain  on  the  diffusion  of 
NH3  is  assumed  to  be  the  same  as  that  for  H2O. 
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Table  II.  Summary  of  the  kinetic  parameters  used  in  this  work. 


Parameter 

XL 

‘^A,(N,inc,ox) 

□ 

n0,(N,inc.ox) 

-il~A,(0.inc,ox) 

p 

r'0.(O.inc.ox) 

■^^A,(0,!nc,Si) 

^O.(O.mc.Si) 

•^^A,(N,mc,Si) 

^O.(N.inc.Si) 

\F 

“il~A,(0,hop.ox) 

p 

0.(0. hop. ox) 
•^A.(N.hop.ox) 

p 

r'0,(N,hop.ox) 


T>900‘C 
2x1  O'3 
0.75 

1.6  eV 

1.2x106  sec-1 

1.4  eV 

1 .2x1  Q6  sec'1 
2.05  eV 
1012  sec'1 
2.05  eV 
1012  sec'1 
0.78  eV 
2.2x1 09  sec'1 
1.2  eV 

3.5x1 07  sec'1 


T<900‘C 
2. 7x1  O'3 
0.75 

1.6  eV 

1 .2x1 06  sec'1 

1.4  eV 

1 .2x1 06  sec'1 

1.6  eV 
1012  sec'1 

1.6  eV 
1 012  sec'1 
1.17  eV 
2.2x1 09  sec'1 

1.4  eV 

3.5x1 07  sec'1 


62 


The  values  of  the  kinetic  parameters  used  in  the  simulations  are 
summarized  in  Table  II. 

IV.  D.  Results  and  Discussion 

Figure  3.12  shows  results  of  the  kinetic  modeling  for  a  1000'C  nitridation 
temperature  at  various  times,  for  the  simplest  case,  in  which  the  nitridation  is 
irreversible  and  there  are  no  strain  effects  (/'.  e.  the  kinetic  parameters  are  the 
same  in  each  dielectric  layer).  The  peak  N  concentration  near  the  interface 
forms  very  early  in  the  nitridation  process  as  the  substrate  reacts  with  the 
diffusing  species,  which  initially  consist  primarily  of  the  diffusing  N  species. 
The  abrupt  jump  in  the  interfacial  nitrogen  concentration  shown  in  Fig.  3.12 
corresponds  to  the  discontinuity  in  the  nitrogen  incorporation  probability  at 
the  SiC^/Si  interface.  As  the  amount  of  nitrogen  incorporated  in  the  oxide 
increases,  more  diffusing  oxygen  becomes  available  with  increasing 
nitridation  time.  The  increasing  concentration  of  the  diffusing  O  species 
relative  to  the  diffusing  N  species  results  in  a  more  oxygen-rich  oxynitride 
forming  at  the  interface  as  the  reaction  proceeds.  In  addition,  the  reaction  of 
the  substrate  with  the  diffusing  species  means  that  the  interface  moves  away 
from  the  initial  interfacial  nitrogen  peak  with  increasing  nitridation  time. 
These  two  effects  result  in  a  peak  in  the  interfacial  nitrogen  concentration 
which  appears  to  move  away  from  the  interface  with  increasing  nitridation 
time,  as  shown  in  Fig.  3.12.  Thus,  even  the  simplest  model  qualitatively 
reproduces  observed13’1 4-17  nitrogen  distributions  for  the  same  nitridation 
conditions.  However,  there  are  quantitative  discrepancies.  The  width  of  the 
calculated  interfacial  nitrogen  peak  is  much  less  than  the  observed  width  (see 
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Fiqure3  12.  Calculated  distribution  of  N  in  a  SiOxN  film  produced  by 
nitriding  a  92  A-thick  SiO.  film  at  1000'C  for  (d)  5  minutes,  (b)  10 
minutes,  (c)  20  minutes,  and  (d)  30  minutes. 
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Fig.  3.3).  The  broadening  of  the  interfacial  peak  in  nitrogen  concentration 
may  be  a  manifestation  of  interfacia!  strain. 

The  effects  of  interfaciai  strain  on  the  calculated  nitrogen  distribution  for  a 
30  minute  mtridation  at  1000‘C  are  shown  in  Fig.  3.13.  The  calculated  curve 
from  Fig.  3.12(d)  (no  strain  effects)  is  shown  in  Fig.  3.13(a).  Figure  3.1 3 f b : 
shows  the  effect  of  including  a  decreased  activation  energy  for  nitrogen 
incorporation  in  the  oxide  near  the  interface,  as  discussed  in  the  previous 
section.  Although  some  broadening  of  the  interfaciai  peak  in  nitrogen 
concentration  is  achieved,  the  discontinuity  due  to  the  dielectric/substrate 
interface  is  still  clearly  evident.  Figure  3.13(c)  shows  the  effect  of  including  a 
decreasing  activation  energy  for  diffusion  near  the  interface,  in  conformity 
with  the  bond  strain  gradient  model.  Including  the  enhanced  diffusivity  near 
the  interface  brings  the  calculated  curve  into  closer  agreement  with  the 
measured14'17  nitrogen  distribution,  and  , educes  the  magnitude  of  the 
discontinuity  in  the  nitrogen  concentration.  A  slight  further  improvement  can 
be  obtained  by  extending  the  region  of  graded  reactivity  to  the  substrate  side 
of  the  interface.  Figure  3.13(d)  shows  the  calculated  nitrogen  distribution 
which  is  obtained  by  reducing  the  probability  (increasing  the  activation 
energy)  for  incorporation  of  nitrogen  and  oxygen  in  the  two  silicon  layers 
closest  to  the  interface.  The  calculated  curves  from  Figs.  3.13(a)  and  3.13(d) 
are  shown  in  Fig.  3.14  together  with  the  measured  nitrogen  distribution  from 
Fig.  3.3(b). 

While  good  agreement  can  be  obtained  between  the  experimental  and 
calculated  nitrogen  distributions  for  a  nitridation  temperature  of  1000‘C,  a 
critical  test  of  the  essential  correctness  of  the  kinetic  model  is  its  ability  to 
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Figure  3.13.  (a)  Calculated  curve  from  Fig.  3.12(d).  (b)  Effect  of  including  an 
activation  energy  for  nitrogen  incorporation  in  the  oxide  which 
decreases  linearly  within  26  A  of  the  interface  as  the  interface  is 
approached,  (c)  Same  as  (b),  but  including  activation  energies  for 
hopping  of  diffusing  species  in  the  oxynitride  which  decrease  linearly 
in  the  mterfacia!  region,  (d)  Same  as  (c)  but  extending  the  region  of 
graded  reactivity  to  include  the  two  Si  layers  in  the  substrate  which  are 
closest  to  the  SiOxNy/Si  interface. 
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Figure  3.14.  Comparison  of  the  experimental  nitrogen  distribution  in  a 
SiOxNy  film,  produced  from  a  100  A-thick  SiO?  nitrided  for  30 
minute's  at  1000'C  (squares,  from  Fig.  3.3(b)),  to  the  calculated 
nitrogen  distributions  (a)  from  Fig.  3.13(a),  and  (b)  from  Fig.  3.13(d). 


predict  nitrogen  distributions  obtained  at  other  nitridation  temperatures. 
Figure  3.15(a)  shows  the  calculated  nitrogen  distribution  for  a  10  minute 
nitridation  at  li50‘C,  using  the  same  kinetic  parameters  that  were  used  to 
calculate  the  nitrogen  distribution  shown  in  Figs.  3.13(d)  and  3.14(b).  The 
experimental  nitrogen  distribution  from  Fig.  3.4(b)  is  also  shown  in  Fig.  3.15. 
Reasonable  agreement  between  the  experimental  data  and  the  simulation  is 
evident.  The  nitridation  occurs  on  a  faster  time  scale  due  to  the  higher 
temoerature,  and  the  peak  interfacial  nitrogen  concentration  saturates  at  a 
lower  value,  as  can  be  seen  by  comparing  Figs.  3.14  and  3.15.  The  lower 
value  of  the  peak  interfacial  nitrogen  concentration  is  due  to  the  increasea 
probability  of  NH3  reacting  with  Si-0  bonds  at  the  elevated  temperature, 
resulting  in  more  nitrogen  incorporating  in  the  oxynitride  and  less  diffusing  to 
the  interface.  The  increased  nitrogen  incorporation  in  the  oxynitride  also 
results  in  an  increase  in  the  concentration  of  diffusing  oxygen,  so  that  a  more 
oxygen-rich  oxynitride  forms  at  the  interface  compared  to  the  1000'C 
nitridation.  The  fact  that  the  simulation  correctly  predicts  this  key  feature 
strongly  suggests  that  the  essential  physics  of  the  nitridation  process  is 
contained  in  the  kinetic  model. 

Figure  3.15(b)  shows  the  calculated  nitrogen  distribution  which  is 
obtained  using  the  same  kinetic  parameters  used  in  the  calculation  of  Fig. 
3.15(a),  but  also  including  the  effects  of  a  reversible  nitridation  reaction,  i.e. 
the  reverse  of  reaction  (3.4)  is  allowed  to  occur,  as  discussed  in  Section  IV. 
Somewhat  better  agreement  between  experiment  and  the  simulation  is 
obtained  assuming  the  nitridation  is  reversible.  However,  this  assumption 
results  in  a  calculated  nitrogen  distribution  at  a  nitridation  temperature  of 
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Figure  3.15.  Comparison  of  the  experimental  nitrogen  distribution  in  a 
SiOxNv  film,  produced  from  a  100  A-thick  SiO?  film  nitrided  for  10 
minutes  at  1150'C  (squares,  from  Fig.  3.4(b)).  to  the  calculated 
nitrogen  distribution  (a)  using  the  same  kinetic  parameters  as  in  Fig. 
3.13(d),  and  (b)  including  the  effects  of  a  reversible  mtridation. 
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1000'C  which  agrees  less  well  with  the  experimental  nitrogen  distribution,  as 
shown  in  Fig.  j.16.  The  data  shown  in  Figs.  3.14  and  3.15  can  be  seen  to  be 
most  consistent  with  an  irreversible  nitridation.  The  results  ot  the  kinetic 
mode!  thus  suggest  that  the  energy  barrier  for  the  reverse  of  reaction  ;2  t, 
has  been  increased  relative  to  the  forward  reaction,  making  oxidation  of  a 
mtrided  oxide  less  likely.  This  suggestion  is  consistent  with  the  known 
oxidation  resistance  of  nitrided  oxides. 50'52'53'60 

The  oossibility  that  the  nitridation  kinetics  are  different  above  anc  below 
'•900‘C  introduces  more  uncertainty  into  the  calculation  of  nitrogen 
distributions  at  low  temperatures.  Figure  3.17(a)  shows  the  calcuiatec 
nitrogen  distribution  for  a  60  minute  nitridation  at  a  nitridation  temperature  of 
800*C,  using  kinetic  parameters  which  are  appropriate  for  low  temperature 
nitridation,  as  discussed  in  Section  IV. C.  Also  shown  in  Fig.  3.17  are  the 
experimental  nitrogen  nitrogen  distributions  from  Fig.  3.5(a)  for  two  fi.ms 
mtrided  at  the  same  conditions.  For  comparison,  Fig.  3.17(b)  shows  the 
calculated  nitrogen  distribution  using  the  same  kinetic  parameters  usee  at 
Higher  temperatures.  Both  sets  of  parameters  result  in  calculated  nitrogen 
distributions  which  predict  a  key  feature  of  the  measured  nitrogen  distribution: 
the  peak  interfacial  nitrogen  concentration  occurs  at  the  interface  and  does 
not  occur  away  from  the  interface  for  long  nitridation  times,  as  at  higher 
temoeratures.  At  this  iow  nitridation  temperature,  the  probability  of  reaction 
!'3.4t  occurring  is  iowered  so  that  more  nitrogen  diffuses  to  the  dielectric' 
substrate  interface.  Sufficient  nitrogen  reacts  with  the  substrate  to  form  a 
affusion  carrier  before  a  significant  amount  of  diffusing  oxygen  is  released  by 
me  slower  substitution  reaction  occurring  in  the  oxide.  The  quantitative 
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gure  3.16.  ComDarison  of  the  experimental  data  points  from  Fig.  3.3(b) 
(squares)  to  the  calculated  nitrogen  distribution  using  the  same  kinetic 
parameters  as  in  Fig,  3.15(b)  for  a  30  minute  mtndation  at  1000'C. 
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Figure  3.17.  Comparison  of  the  experimental  nitrpgen  distribution  in  SiO  N 
films,  produced  from  100  A-thick  Si02  films  nitrided  for  60  minute'"'  at 
800’C  (solid  and  open  circles,  from  Fig.  3.5(a)),  to  the  calcuib.  I 
nitrogen  distribution  (a)  using  kinetic  parameters  appropriate  for  low 
temperature  mtridation,  and  (b)  using  kinetic  parameters  appropriate 
for  high  temperature  mtridation. 


ciscrepancy  between  Fig.  3.17(b)  and  the  measured  nitrogen  distriouticn 
suggests  that  the  mtridation  kinetics  are  indeed  different  above  and  beiow 
~F0Q‘C.  However,  the  quantitative  discrepancy  between  rig.  o. 17(a)  and  the 
measured  nitrogen  distribution  is  an  indication  that  the  low  temperature 
Kinetic  parameters  are  not  optimized.  The  most  iikeiy  source  of  error  is  the 
assumption  that  R0  for  the  various  reactions  does  not  vary  with  temperature, 
even  though  AEA  does  vary. 

The  calculated  nitrogen  distributions  shown  in  Figs.  3.12  to  3.17  were  a! 
using  an  initial  Si02  thickness  near  100  A.  The  effect  of  increasing  oxice 
thickness  is  shown  in  Fig.  3.18  for  a  30  minute  mtridation  at  1000'C.  The 
peak  interfacial  nitrogen  concentration  is  seen  to  decrease  in  magnitude  anc 
occur  closer  to  the  interlace  with  increasing  initial  oxide  thickness.  This  trend, 
together  with  the  reduced  interfacial  nitrogen  incorporation  at  higner 
temoeratures  previously  discussed,  may  explain  the  failure  to  observe50  an 
interfacial  nitrogen  buildup  for  a  500  A  oxide  mtrided  at  1200"C.  The 
magnitude  cf  the  decrease  shown  in  Fig.  4.18  is  likely  to  be  larger  than  what 
would  occur  experimentally,  howe''er,  since  an  interfacia'  nitrogen  buildup  is 
clearly  evident67'69  in  400  A  thick  oxides  which  have  been  nitrided  at  1 100’C. 

In  addition  to  the  quantitative  discrepancies  noted  above  for  low 
.mtridation  temperatures,  quantitative  discrepancies  between  the  data  and  the 
results  of  the  simulations  are  apparent  in  the  higher  temperature  mtridat'ons 
snown  m  Figs.  3.14  and  3.15.  The  nitrogen  concentration  m  the  bulk 
: between  the  surface  and  interfacial  peaks  in  nitrogen  concentration)  of  the 
‘:;m  's  generally  underestimated,  especially  for  long  mtridation  times.  In  *he 
ovulations,  a*ter  sufficient  nitrogen  has  incorporated  m  the  surface  ’aver  of 
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Figure  3.18.  Calculat°  ,ogen  distribution  using  the  same  kinetic 
parameters  a~  .  u.  13(d)  for  a  30  minute  nitridation  at  1000'C  for  an 
initial  oxide  thickness  of  (a)  100  A,  (b)  130  A,  and  (c)  200  A. 


the  mtrlded  oxide,  no  additional  diffusing  nitrogen  can  enter  the  oxynitride  film 
from  the  gas  phase  and  the  distribution  of  nitrogen  in  the  film  reaches  the 
steady  state.  Experimentally,  however,  nitrogen  incorporation  continues  in 
the  bulk  of  the  film  even  after  the  N  distribution  in  the  interfacial  region  and 
the  surface  N  concentration  cease  to  change. 1“'17’67'69  This  may  indicate  that 
the  heavily  nitrided  surface  layer  is  too  thin  to  be  an  effective  diffusion  oarrier, 
whereas  in  the  simulations  oniy  a  single  layer  (2.63  A)  of  forms  a 

perfect  diffusion  barrier.  Another  discrepancy  is  apparent  in  the  interfacial 
region,  where  the  simulations  show  the  nitrogen  concentration  decreasing 
more  rapidly  than  is  experimentally  observed  in  the  region  between  the 
interface  and  the  maximum  in  the  interfacial  nitrogen  concentration. 
Preferential  etching  could  account  for  this  discrepancy  if  the  distribution  of 
nitrogen  is  laterally  nonuniform,  since  the  available  experimental  evidence 
shows  that  a  higher  nitrogen  concentration  lowers  the  etch  rate.  Nitrogen-rich 
regions  would  then  etch  last,  skewing  the  distribution  towards  the  interface. 
Interfacial  kinetics  which  differ  from  that  assumed  in  this  simple  model  may 
also  cause  the  discrepancy. 

V.  Spectroscopic  Eliipsometry  Measurements 

Spectroscopic  eliipsometry  measurements  were  taken  on  the  samples 
which  were  prepared  as  summarized  in  Table  I.  Two  different  single-layer 
models  were  fit  to  data  taken  on  sample  2.  in  the  first,  the  SiOxNy  layer  was 
modeled  as  a  mixture  of  Si02  and  voids.  The  best  fit  obtained  had  a  mean 
sauare  error  (MSE)  of  5.C  „  1 0"4  with  an  oxide  thickness  tox  =  127  A  when  the 
void  fraction  fy  was  fixed  at  0.  When  fv  was  allowed  to  vary,  the  best  fit  (MSE  = 


3.6  x  1 0'4)  was  obtained  with  tox  =  1 1 2  A  and  fv  =  -0.35. 

The  second  model  considered  the  SiOxNy  layer  to  be  a  physical  mixture 
of  Si02  anc  SigN,^.  It  is  to  be  noted  that  in  this  special  case,  a  pnysical 
mixture  is  eouivalent  to  a  chemical  mixture.  This  is  because  for  both  S1O5 
and  S13N4  absorption  in  the  energy  range  of  interest  is  dommatea  by  the 
ione  pair  electrons  on  0  or  N  48  rather  than  the  bonding  electrons.  Within  the 
second  model,  the  best  fit  (MSE  =  1.3x1  O'4)  was  obtained  for  tox  =  1 1 6  A  and 
a  nitride  fraction  fN  =  0.065.  For  comparison,  the  MSE  between  two  sets  of 
data  taken  on  the  same  sample  was  1.4x1  O'4. 

Samples  1  to  9  can  clearly  provide  no  information  on  the  interface,  since 
the  data  can  be  well  fit  with  a  simple  one-layer  model.  The  lack  of  sensitivity 
can  be  traced  to  the  thin  dielectric  film.  For  a  thickness  of  100  A,  the  phase 
difference  between  light  reflected  at  the  ambient/film  and  film/substrate 
interfaces  is  insufficient  to  cause  the  extinction  of  s-polarized  light  which  gives 
rise  to  a  peak  in  tan  y,  and  hence  enhanced  interfacial  sensitivity. 
Nevertheless,  a  two-layer  model  has  been  used64  on  spectroscopic 
ellipsometry  data  taken  on  200-A-thick  nitrided  oxides.  However,  fv  was  the 
only  new  information  measured,  with  the  N  fraction  provided  by  RBS.  For 
these  thin  dielectric  films,  spectroscopic  ellipsometry  is  clearly  insufficiently 
sensitive  to  determine  the  complex  nitrogen  distributions  observed  in  nitrided 
oxides. 

Figures  3.19  and  3.20  show  the  data  taken  on  samples  10  and  1 1 ,  as 
well  as  a  samDie  taken  from  a  wafer  oxidized  at  the  same  time  under  the 
same  conditions,  but  with  no  nitridation.  Comparing  the  oxide  curve  in  Figs. 
3.19(a)  and  3.20(a)  to  similar  data  presented  elsewhere,47  the  major 


qualitative  discrepancy  is  that  in  Fig.  3.20(a),  the  minimum  in  cos  A  near  2.6 
eV  does  not  reach  -1.  Such  an  observation  has  previously1 03  been 
interpreted  as  eviaence  of  thickness  nonuniformity.  The  degree  of  thickness 
vacation  required,  however,  is  several  tens  of  Angstroms.  This  explanation 
seems  unlikely  for  the  device-auality  thermal  oxide  sample  used.  As  seen  m 
Figs.  3.20(b)  and  (c),  this  situation  becomes  worse  with  an  increasing  aegree 
of  nitridation.  Two  possible  causes  are  roughening  of  either  the  surface  or 
interface,  or  laterally  nonuniform  nitnonen  incorporation  in  the  oxide.  In  this 
latter  case,  nitrogen-rich  regions  would  have  a  higher  index  of  refraction  and 
would  appear  to  be  optically  thicker  than  surrounding  regions. 

Several  models  were  fit  to  the  data  taken  on  sample  11,  shown  in  Figs. 
3.19(c)  and  3.20(c).  A  good  fit  could  be  obtained  to  tan  y  in  all  cases,  so  only 
the  fit  to  cos  A  wiil  be  shown.  Figure  3.21  shows  the  result  of  fitting  a 
single-layer  model  to  the  data,  consisting  of  a  physical  mixture  of  Si02  and 
Si3N4.  The  best  fit  (MSE  =  6  x  10'3)  was  obtained  for  tox  =  1074  A  and  fN  = 
0.034.  As  can  be  seen,  cos  A  is  not  fit  well  near  2.5  eV.  A  fit  to  cos  A  only 
yields  MSE  =  6.65  x  10'4,  tQX  =  1033  A,  and  fN  =  0.091.  However,  in  this  case 
tan  y  has  a  maximum  value  of  19.5,  compared  to  the  observed  value  near  12 
in  the  data.  Including  additional  layers  in  the  model  at  the  surface  and  near 
the  interface  with  each  layer  allowed  to  have  different  values  of  fN  only 
marginally  decreased  the  MSE  (by  <  10%). 

Figures  3.22  and  3.23  shows  data  taken  on  the  same  sample  for  two 
different  aoerture  sizes.  This  technique  was  previously  used103  to 
demonstrate  the  effects  of  thickness  nonuniformities.  The  aperture  sizes  used 
r'0r  the  aperture  nearest  the  detector)  were  0.6  mm  and  2  mm.  This 
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Figure  3.21 .  Best  fit  to  the  data  from  Figs.  3.19(c)  and  3.20(c),  obtained  using 
a  single-layer  model  in  which  the  dielectric  consists  of  a  physical 
mixture  of  SiCU  and  Si3N4.  Only  the  fit  to  cos  A  is  shown,  since  tan  v 
could  be  fit  wen  enough  that  no  difference  between  the  data  and  the  fit 
is  apparent  to  the  eye. 


cos  A 


corresponds  to  a  smaller  spot  on  the  sample,  since  the  beam  diverges 
between  the  sample  and  the  detector.  Some  attempts  were  maoe  to  include 
the  effects  of  nonuniformities  in  the  modeling.  Within  a  single-layer  model,  a 
Gaussian  distribution  of  thicknesses  about  some  mean  value  o'  t  was 

C  x 

assumed.  The  oest  fit  was  for  a  Gaussian  full  width  at  half  maximum  of  26  A. 
out  the  quality  of  the  fit  was  not  better  than  that  shown  in  Fig.  3.21 .  A  similar  fit 
was  attempted  assuming  nonuniformities  in  fN,  but  no  improvement  in  the  fit 
could  be  obtained.  It  may  be  that  the  nonuniformities  are  greater  than 
assumed,  and  the  fits  are  not  converging  at  the  true  minimum  in  the  MSE.  At 
this  time,  the  data  is  not  understood. 

Seme  model  data  was  calculated  in  order  to  determine  limits  on  the 
usefulness  of  spectroscopic  ellipsometry  in  the  study  of  nitrided  oxides,  and  to 
determine  the  usefulness  of  infrared  ellipsometry  in  studying  this  material 
system.  Figures  3.24  and  3.25  shows  calculated  spectra  for  a  1030  A  thick 
film  on  Si  with  fN  -  0.1  in  the  dielectric.  As  can  be  seen,  tan  v  is  a  slowly 
varying  function  below  1.6  eV,  and  cos  A  is  rapidly  varying.  The  MSE 
between  these  curves  and  similar  curves  calculated  with  a  15-A-thick  SigN^ 
layer  at  the  interface  and  reducing  the  bulk  dielectric  layer  thickness  by  the 
same  amount  is  0.15.  However,  the  differences  are  apparent  primarily  in  the 
visible-UV  region.  There  would  thus  appear  to  be  no  advantage  to  using 
infrared  elliDSometry  on  this  material  system 

The  calculated  curve  with  the  15-A-thick  SigN4  interface  layer  and  a 
101 5-A-thick  bulk  with  fN  =  0.1  was  fit  with  a  single-layer  model.  The  best  fit 
had  a  MSE  =  5  x  1 0'4 ,  so  the  interface  is  detectable.  The  two-layer  calculated 
cu rve  was  ccmcared  to  a  ;hree-layer  calculated  curve  in  which  the  interface 
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F;gure  3.24  Calculated  spectrum  of  tan  v  for  a  1030  A-thick  dielectric  on  Si, 
ccnsistma  of  a  physical  mixture  of  10%  Si3N4  and  90%  Si02,  in  the 
energy  range  0.5-5. 7  eV. 
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Figure  3.25.  Calculated  spectrum  of  cos  A  for  a  1030  A-thick  dielectric  on  Si 
consisting  of  a  physical  mixture  of  10%  Si3N4  and  90%  Si02,  in  th< 
energy  range  0.5-5. 7  eV. 


consisted  of  15  A  of  S13N4  on  top  of  15  A  of  S1O2.  The  MSE  was  3.7  x  10‘3 
However,  when  a  two-layer  mode!  was  fit  to  the  three-layer  calculated  curve 
the  MSE  was  5  x  10'6.  well  below  the  noise  level.  Spectroscopic  elliDSomet7 
would  thus  not  be  capable  of  distinguishing  a  nitripe-rich,  interface  from  a 
nitrogen-rich  region  separated  from  the  interface  by  an  oxygen-rich  region. 
Beth  of  these  structures  occur  in  the  mtrided  oxides,  depending  on  tne 
mtridation  conditions.  Spectroscopic  ellipsometry  wouid  thus  aopear  tc  be  of 
limited  usefulness  in  studying  nitrided  oxides. 

VI.  Surrmary  and  rtnnr!nsions 

In  summary,  the  thermal  nitridation  of  SiC>2  films  on  Si  has  been  stuciec 
by  measuring  the  nitrogen  concentration  depth  profiles  with  XPS  as  a 
function  of  oxidation  temperature  and  nitridation  time  and  temperature.  The 
peak  interfacial  nitrogen  concentration  initially  occurs  at  the  S i 0 x N y / S i 
interface,  then  appears  to  move  away  from  the  interface  with  increasing 
nitridation  time  for  sufficiently  high  niiiiciation  temperatures.  The  behavior  of 
AVfb|  correlates  well  with  the  observed  nitrogen  distributions.  This 
correlation  is  shown  to  consistently  follow  from  the  simple  postulate  that  it  is 
the  strain-dependence  of  the  kinetics  of  formation  of  defect  states  which  allow 
lowering  of  the  system  energy  by  relieving  interfacial  strain.  The  distribution 
of  defects  is  deduced  from  the  intensity  of  a  fluorine  marker,  and  is  found  to 
be  consistent  with  the  proposed  model  explaining  the  behavior  of  flatband 
voltage  shifts  as  a  function  of  nitridation  conditions.  Macroscopic  intrinsic 
stress  is  not  observed  to  affect  the  nitrogen  distribution,  but  does  affect  the 
defect  formation  kinetics. 
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Modeling  of  the  kinetics  of  the  nitndatio n  of  S1O2-S1  in  NH3  has  oeen 
done  in  this  study,  and  reveals  a  mechanism  in  which  diffusing  species, 
initially  consisting  primarily  of  the  nitrogen  species,  react  with  the  substrate. 
formmg  a  nitrogen-rich  oxynitride  a’  the  interface.  At  sufficiently  h:a~ 
mtrication  temperaturt-s.  an  oxygen-rich  oxynitride  subseauently  forms  at  the 
interface  due  to  reaction  of  the  substrate  with  an  increasing  concentration  of 
diffusing  oxygen  which  has  been  displaced  by  the  slower  mtridation  of  the 
dielectric.  The  data  are  consistent  with  the  mechanism  discussec  above, 
orovided  the  influence  of  the  mterfacia!  strain  on  the  nitncation  ana  oxicaticn 
kinetics  ■?  evpii-rt!,  amounted  for.  The  fact  that  the  XPS  data  Qualitatively 
confirm  the  predictions  of  the  kinetic  model  is  an  indication  that  the  essential 
physics  of  the  nitridation  process  is  contained  in  the  model,  although 
quantitative  discrepancies  indicate  that  the  kinetic  parameters  are  not  yet 
optimized.  In  spite  of  the  discrepancies,  the  model  successfully  predicts  that 
(a)  the  maximum  nitrogen  concentration  in  the  interfacial  region  initially 
occurs  at  the  interface  and  occurs  away  from  the  interface  with  increasing 
mtridation  time  at  sufficiently  high  mtridation  temperatures,  (b)  the  maximum 
nitrogen  concentration  in  the  interfacial  region  remains  at  the  interface  even 
for  long  mtridation  times  at  low  nitridation  temperatures,  and  (c)  the  peak 
interfacial  nitrogen  concentration  is  lower  for  the  highest  mtridation 
temperature  studied.  The  simulations  provide  valuable  insight  into  the 
reasons  for  these  exoerimentally  confirmed  observations,  and  also  provice 
insight  into  the  role  of  interfacial  strain  in  the  nitridation  kinetics. 
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Chapter  4 

Studies  of  Semiconductcr-on-lnsulator  (SOI)  Structures 


I.  Introduction 

Semiconductor-on-insulator  (SO!)  structures  consisting  of  Si/SiOo.Si  are 
a  potentially  low-cost  alternative  to  silicon-on-saphire  structures.  Such 
structures  show  promise  in  high-speed  complementary  metal-oxide-semicon- 
ductor  (CMOS',  devices  with  applications  in  radiation  hard  electronics  an c 
VLSI  circuits.  The  formation  of  SCI  structures  can  be  accomplished  in  several 
ways.  Two  techniques  have  been  studied  as  part  of  this  task:  ion 
implantation  and  zone-melting  recrystallization  (ZMR). 

The  formation  of  SOI  structures  by  ion  implantation  invo.ves  the  hig- 
dosage  implantation  of  ions  into  a  single  crystal  Si  substrate.  Heating  the 
substrate  during  the  implantation  (typically  to  500’  C)  aids  both  the 
redistribution  of  the  implanted  species  and  the  annealing  of  implantation 
damage  in  the  Si  surface  layer  and  the  Si  substrate.  A  subsequent  high 
temDerature  anneal,  during  which  the  implanted  species  react  with  the  Si, 
forms  the  buried  insulator  layer.  The  implanted  species  can  be  oxygen  to 
form  a  buried  Si02  layer, 19'22-106'108  nitrogen  to  form  a  buried  S^^ 
layer,109'112  or  both  species  to  form  a  buried  SiOxNy  layer.113  The  post- 
implantation  annealing  temperature  is  typic  ’  1100’  C  to  1300’  C.  though 

recently  excellent  results  have  been  obtained  at  1405’  C.107  In  addition  to 
chemical,  electrical,  and  morphological  characterizations,19'22'106707'109'" 13 
tr.s  formation  of  SOI  structures  by  ion  implantation  has  also  been  modeled.108 


in  this  work,  the  morphology  of  the  buried  SiC^/Si  interfaces  has  been 
stuoied  by  SEM,  and  the  chemical  structure  of  the  interface  of  thermally 
Grown  S1O9  on  the  surface  Si  has  been  studied  by  XPS.  These  studies  are 
described  in  section  II. 

The  formation  of  SOI  structures  by  ZMR  involves  the  deposition  of 
pciycrystalli ne  Si  onto  a  S1O2  layer  which  has  been  thermally  grown  on  3 
single  crystal  Si  substrate.  A  S1O2  capping  layer  is  deposited,  and  u  graphite 
strip  heater  is  passed  over  the  structure  to  melt  the  poly-Si  layer,  which  wnen 
cooied  recrystallizes  to  form  a  single  crystal  Si  layer.  Experimental 
measurements23  and  modeling114  of  the  formation  of  such  structures  have 
been  done.  It  has  been  reported23  that  nitridation  of  the  capping  layer 
improves  the  quality  of  the  recrystailized  layer.  It  has  been  suggested23  that 
the  introduction  of  nitrogen  at  the  capping  SiC^/deposited  Si  interface 
improves  the  wetting  of  the  capping  iayer  by  the  Si  when  it  is  melted. 
However,  attempts  to  detect  nitrogen  at  the  capping  SiC^/depcsited  Si 
interface  have  not  been  successful.23  In  this  work,  XPS  is  used  in  an  attempt 
to  detect  the  possible  presence  of  nitrogen  at  the  interface.  The  results  of 
these  studies  is  reported  in  section  III. 

II.  SOI  Formed  by  Ion  Implantation 

Four  SOI  samples  formed  by  ion  implantation  were  obtained  from  Dr. 
Peter  Hemrrent  of  the  University  of  Surrey,  England.  The  sample  preparation 
is  described  elsewhere.19'22  All  of  the  samples  were  implanted  with  a 
desage  of  1.8  x  1018  cm-2  oxygen  ions  with  an  implantation  energy  of  200 
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keV.  The  implantation  temperature  for  two  of  the  samples  was  500‘C.  One 
was  suDsequently  annealed  at  1150’  C  for  two  hours,  while  the  other 
recieved  no  further  processing.  The  purpose  of  studying  these  two  samples 
is  to  determine  the  effects  of  annealing.  The  other  two  samples  were  botn 
annealed  at  1 1 50*  C  for  two  hours;  the  implantation  temperature  for  one  was 
425'  C  ana  for  the  other  was  525'  C.  The  purpose  of  studying  these  two 
samples  is  to  determine  the  effect  of  varying  implantation  temperature. 

Prior  to  any  further  processing,  the  surfaces  of  ail  four  samples  were 
studied  using  XPS  and  SEM.  Both  the  XPS  and  SEM  data  suggest  that  the 
surface  Si  film  is  good  quality  single  crystal  Si,  as  discussed  below.  The 
SEM  micrographs  show  no  evidence  of  grains,  which  would  have  been  an 
indication  of  polycrystalline  material,  suggesting  that  the  surface  Si  layer  is 
either  single  crystal  or  amorphous.  However,  SEM  micrographs  taken  at 
glancing  angles  (70'  from  the  normal  to  the  sample  surface)  do  show 
evidence  of  varying  degrees  of  surface  roughness,  apparent  as  small 
"hillocks"  on  the  surface  approximately  100  A  in  diameter,  as  shown  in  Fig. 
4.1.  The  XPS  data  show  substantial  differences  in  the  native  oxides  on  these 
samples.  Figure  4.2  shows  the  Si  2 p  region  of  the  photoelectron  spectrum, 
showing  both  the  substrate  (near  100  eV)  and  oxide  (near  104  eV)  peaks. 
The  substrate  peak  shape  is  consistent  with  the  surface  Si  layer  being  single 
crystal,  and  does  not  show  the  broadening  exhibited44  by  amorphized  Si. 
Differences  in  oxide  thicknesses  are  exhibited  in  the  differing  oxide/substrate 
intensity  ratios,  and  differences  in  native  oxide  stoichiometry  are  exhibited  in 
the  differing  oxide-to-substrate  energy  separations  and  oxide  oeak  shapes.  It 
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Figure  4.1.  SEM  microqraphs  of  the  surfaces  of  SOI  samples  formed  by  ion 
implantation  of  oxygen  at  (a)  500°C,  no  anneal,  (b)  425°C,  2  hour 
1150°C  post-impiantation  anneal,  (c)  500°C,  2  hour  1150°C  anneal, 
and  (d)  525°C.  2  hour  1 150°C  anneal. 
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Figure  4.2.  XPS  spectra  ot  the  Si  2p  region  for  the  samples  corresponding  to 
'he  data  in  Fig.  4.1. 
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has  been  suggested1 15  that  the  suboxide  distribution  at  the  SiC^/Si  mtedace 
is  affected  by  the  surface  step  density.  The  native  oxice  stoichiometry  may 
thus  prove  useful  as  a  measure  of  the  perfection  of  the  crystalline  surface. 

The  buriea  S1O2  layer  has  also  been  studiea  with  both  XPS  ana  SEM. 
The  buried  oxide/surface  Si  and  buried  oxide/substrate  interfaces  were 
exposed  by  etching  away  the  surface  Si  layer  and  buried  oxide  layer 
respectively.  The  etching  technique  used  for  removing  the  surface  Si  layer 
was  the  XeF2  etching  described  briefly  in  Chapters  2  and  3  and  in  mce 
detail  elsewhere.37-38  The  SEM  micrographs  of  the  buried  SiO?- surface  S 
interface  exposed  in  this  way  are  shown  in  Fig.  4.3.  Figure  4.4  shows  the 
corresponding  SEM  micrographs  at  a  higher  resolution.  The  SEM 
micrographs  for  the  unannealed  sample,  shown  in  Figs.  4.3(a)  and  4.4(ai. 

ft 

show  a  smooth  surface,  as  expected  if  the  oxygen  implantation  profile  is 
laterally  uniform.  Figures  4.3(b)  through  (d)  and  4.4(b)  through  (d)  are  the 
SEM  micrographs  showing  the  effect  of  increasing  implantation  temperature. 
Figures  4.3(b)  and  4.4(b)  show  the  roughest  interface,  indicating  that  for  the 
low  implantation  temperature  used  (425'  C),  the  redistribution  of  implanted 
oxygen  is  too  slow  to  form  a  well-defined  buried  oxide  with  sharp  interfaces. 
Figures  4.3  (c)  and  (d)  and  4.4(c)  and  (d)  show  progressively  smoother 
interfaces  with  increasing  implantation  temperature.  The  difference  between 
micrographs  (c)  and  (d)  in  both  Figs.  4.3  and  4.4  seems  surprising  at  first 
giance,  since  the  implantation  temperatures  differ  by  only  25'  C.  However, 
the  implantation  temperature  range  500'  to  550'  C  is  critical  in  determining 
the  nature  of  the  surface  Si/buried  Si02  interface.  Below  500'  C.  the  Si 
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Figure  4.3.  SEM  micrographs  (magnification  10-1 3kx)  of  the  surfaces  of  SCI 
samples  formed  by  ion  imolantation  of  oxygen  at  (a)  500  C.  no  anneal. 
' b)  425  C,  2  hour  11503C  post-implantation  anneal,  (ci  500'C.  2  hour 
1  f  50 "C  anneal,  and  (d)  525CC,  2  hour  1 150°C  anneal.  The  surface  Si 
layer  iias  been  removed  by  XeF2  etching,  exposing  the  buriec  oxide. 


F'Qjre  4  4  Same  as  F-g.  4  3.  out  with  a  magnification  of  65kx 


nearest  that  interface  is  amcrpnous,  wniie  above  5c0‘  C  it  is  best  aescrioec 
as  ae'ec.ed  crystalline.'15  In  view  of  this,  the  differences  apparent  m  Figs.  4.3 
and  4.4  are  net  unreasonable.  The  XPS  data  taken  on  these  samples  were 
net  consicered  reliable  due  to  sample  charging,  a  consequence  of  the 
thickness  of  the  insulator  layer. 

After  the  studies  discussed  above,  the  buried  SiC^/substrate  interface 
was  exposed  by  removing  the  buried  S1O2  layer  with  HF.  The  SEM 
micrographs  shown  in  Figs  4.5  and  4.6  show  the  same  trends  toted  above 
‘or  the  surface  Si/buried  SiC>2  interface,  /.  e.  the  interfaces  are  smoother  at 
higher  implantation  temperatures.  The  XPS  data  show  the  unannea.ee 
sample  to  nave  a  Si  2 p  lineshape  which  is  characteristic  of  amorphous  Si.““ 
while  the  three  annealed  samples  all  exhibited  lineshapes  characteristic  of 
crystalline  Si.  Only  trace  amounts  of  oxygen  (corresponding  to  a  tew  oercent 
of  a  monolayer)  were  detected,  probably  from  residual  surface  contaminants. 
However,  the  samples  appeared  visually  to  have  a  thick  dielectric  layer, 
exhibiting  the  colors  typical  of  optical  interference  effects.  The  optical 
interference  indicates  a  surface  layer  with  an  index  of  refraction  which  differs 
from  that  of  the  underlying  Si  substrate.  Subsequent  XeF2  etching,  however, 
removed  the  remainder  of  the  samples,  indicating  no  continuous  SiC>2  layer 
was  present.  Some  possible  causes  f  the  optically  distinct  surface  layer 
include  surface  rough, less,  a  region  of  defected  Si,  or  a  region  with  a 
relatively  high  concentration  of  oxygen  or  Si02  inclusions. 

In  Devices  based  on  SOI  structures,  it  is  the  quality  of  ihe  thermally  grown 
SiOg-surface  Si  interface  which  is  important  in  determining  aevice 
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4  5  SFM  micrographs  (magnification  50kx)  of  the  surfaces  of  SOI 
cameras  formed  Py  ion  implantation  of  oxygen  at  (a)  500  C.  no  anneal 
;oi  425  C.  2  nour  1150  C  post-implantation  anneal,  ic)  500'C.  2  ncur 
1  *50  C  anneal,  and  id)  525'C,  2  hour  1 1 50cC  anneal.  The  surface  Si 
a yer  nas  been  removed  by  XeF2  etching,  and  the  bur'ed  oxide  aye'- 
nas  been  removed  with  HF,  exposing  the  buried  oxide  .substrate 
me  dace 


characteristics.  Consequently,  the  samples  which  were  implanted  at  425’  C 
anc  525'  C  and  annealed  for  two  hours  at  1150”  C  were  oxidized  in  dry  O2  tc 
give  a  SiC2  layer  approximately  80  A  thick.  This  oxide  was  then  stripped  with 
HP,  ana  an  80-A-thick  SiC>2  layer  was  again  thermally  grown,  whicn  is  the 
standard  processing  for  device  quality  oxides.  These  thermally  grown  oxides 
were  depth  profiled  using  the  sequential  chemical  etching  technicue 
described  in  Chapter  3  and  elsewhere,1'35  with  XPS  spectra  being  recordec 
after  each  etch  step.  The  thermally  grown  SiC^/surface  Si  interfaces  for  the 
two  samples  appear  to  be  identical  within  experimental  error,  which  seems  tc 
be  inconsistent  with  the  previously  noted  differences  in  the  native  oxides. 
However,  it  is  well  known  that  oxidation  tends  to  smooth  the  interface.  These 
results  suggest  that  even  though  the  buried  SiC>2  has  interfaces  with  the  Si 
substrate  and  surface  Si  layer  with  structures  which  depend  fairly  strongly  on 
implantation  temperature,  the  quality  of  the  thermally  grown  Si02  and  the 
Si02/surface  Si  interface  are  not  significantly  affected. 

III.  SOI  Formed  by  Zone-Melting  Recrystallization  (ZMR) 

Two  SOI  samples  formed  by  ZMR  were  obtained  from  Dr.  Chenson  Chen 
of  MIT  Lincoln  Laboratory.  The  sample  preparation  is  described  elsewhere.23 
A  S1O2  layer  2  pm  thick  was  grown  by  steam  oxidation  on  a  single  crystal  Si 
substrate.  A  polycrystalline  Si  layer  0.5  pm  thick  was  deposited,  followed  by 
the  deposition  of  2  pm  of  Si02  as  a  capping  layer.  One  sample  was  used  as 
a  control,  the  other  was  nitrided  in  NH3  at  atmospheric  pressure  at  1100”  C 
for  two  hours,  oxidized  in  dry  O2  for  20  minutes,  and  nitrided  for  an  additional 
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two  hours.  The  Durpose  of  the  intermediate  oxidation  is  to  reoxidize  the 
surface  diffusion  barrier  which  forms  as  a  result  of  the  nitridation  (see  Chapter 
3)  so  that  additional  NH3  is  able  to  enter  the  film  with  additional  nitridation 
time.  It  has  been  suggested23  that  this  procedure  results  in  the  introduction  ot 
N  at  the  capping  Si02/deposited  Si  interface,  which  improves  the  wetting  c* 
the  capping  layer  by  the  molten  Si  in  ZMR  and  results  in  superior 
recrystallized  films.  However,  attempts  to  detect  the  N  at  the  interface  have 
been  unsuccessful.23  The  purpose  of  the  studies  reported  here  is  to 
determine  if  N  has  been  incorporated  at  the  interface  using  XPS,  which  may 
be  more  sensitive  than  the  techniques  previously  employed.23 

In  order  to  confirm  the  hypothesis  that  N  is  incorporated  at  the  capping 
Si02/deposited  Si  interface,  the  interface  must  be  exposed  for  study  by  XPS. 
It  is  unlikely  that  an  HF  etch  through  2  (im  of  Si02  would  be  uniform  enough 
to  come  within  a  few  tens  of  A  of  the  interface  (within  one  photoelectron 
escape  depth)  without  nonuniformities  developing,  making  data  analysis 
difficult.  Exposure  of  the  interface  from  the  Si  side  is  therefore  desirable  in 
this  case.  In  order  to  expose  the  interface,  the  Si  substrate,  thermally  grown 
Si02  layer,  and  recrystallized  Si  layer  must  be  removed.  This  can  be  done 
by  a  XeF2  etch  to  remove  the  Si  substrate,  followed  by  a  HF  etch  to  remove 
the  thermally  grown  Si02  layer,  followed  by  another  XeF2  etch  to  remove  the 
recrystallized  Si  layer.  One  of  Al  was  deposited  on  the  front  of  both  the 
mtrided  and  control  samples  for  structural  support  of  the  capping  layer  after 
removal  of  the  rest  of  the  structure,  and  both  samples  were  bonded  to  a  Au 
platen  using  high-purity  In. 
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It  was  found  that  the  back  c4  the  samples  had  a  thick  layer  of  S1O2 
apparently  originating  during  the  thermal  oxidation  or  the  deposition  of  tne 
caoDing  layer.  This  was  removed  with  HF  pnor  to  exposure  of  the  samples  to 
Xep2-  After  removal  of  the  Si  substrate,  it  was  found  that  the  HF  etcn  hac 
aDparently  also  removed  the  A!  film  (and  probably  the  capping  layer  with  it 
The  failure  of  the  Al  layer  and  In  bonding  to  protect  the  sample  surface 
requires  that  a  method  of  protecting  the  surface  be  found.  This  problem  was 
solved  by  removing  the  backside  SiC>2  layer  only  from  the  central  region  of 
the  sample.  The  XeF2  etching  then  would  leave  a  ring  of  Si  around  the  edge 
of  the  sample,  which  would  contain  the  HF  etchant  during  removal  of  the 
buried  Si02  layer.  Removal  of  the  recrystallized  Si  layer  would  then  expose 
the  capping  Si02.  In  practice,  the  intended  ring  of  Si  was  not  continuous, 
and  sufficient  etchant  "escaped"  the  ring  to  damage  the  protective  Al  layer. 
As  of  this  writing,  the  problem  of  exposing  the  capping  SiC^/recrystallized  Si 
interface  has  not  been  solved. 
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Chapter  5 

Conclusions 

The  techniques  of  XPS,  SE,  and  SEM  have  been  used  to  study  structural 
and  chemical  inhomogeneities  in  several  electronic  materials  and  device 
structures  of  relevance  to  radiation  hard  electronics.  In  many  ways,  XPS  is 
complementary  to  SE  and  SEM,  and  this  fact  has  been  exploited  in  this  work 
to  elucidate  information  on  chemical  and  structural  properties  of  the  materiais 
and  structures  studied.  The  systems  studied  were  MNOS  structures,  silicon 
oxynitride  formed  by  the  thermal  nitridation  of  SiC>2/Si,  and  SOI  structures. 
Studies  of  MNOS  structures  suggest  that  the  effect  of  H2  annealing  is  to  make 
the  SigN^SiC^  interface  less  abrupt  by  causing  interdiffusion  of  silanoi  ana 
silamine  groups  with  subsequent  oxynitride  formation.  Another  effect  of  the 
annealing  appears  to  be  to  relieve  the  strain  at  the  Si02/Si  interface. 

The  thermal  nitridation  in  NH3  of  Si02  films  on  Si  results  in  the  formation 
of  an  oxynitride  film  which  is  inhomogeneous  both  chemically  and 
structurally.  In  addition  to  the  structural  inhomogeneities  inherent  in  an 
amorphous  material,  the  oxide  near  the  Si02/Si  interface  (within  -30  A)  is 
also  known  to  be  structurally  strained  and  distinct  from  the  bulk  SiC>2.  The 
effect  of  this  strain  on  the  nitridation  kinetics  has  been  shown  to  affect  the 
observed  nitrogen  distribution  in  the  SiOxNy  film.  The  nitrogen  concentration 
depth  profiles  as  a  function  of  oxidation  temperature  and  nitridation  time  and 
temperature  have  been  measured  using  XPS  intensity  analysis  in 
conjunction  with  chemical  etching.  For  this  material  system,  SE  was  found  to 
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be  of  limited  usefulness,  providing  only  some  evidence  of  interface 
roughening  as  a  result  of  nitridation.  However,  this  effect  could  not  be 
quantified  in  this  work.  The  XPS  data  show  that,  at  a  nitridation  temperature 
of  1000'  C,  the  peak  interfacial  nitrogen  concentration  initially  occurs  at  the 
SiOxNy/Si  Interface.  then  appears  to  move  ~20  A  away  from  the  interface  with 
nitridation  Ume  increasing  from  10  to  30  minutes.  For  nitridation  times  beyond 
30  minutes,  the  distribution  of  nitrogen  within  ~30  A  of  the  interface  does  not 
change,  but  nitrogen  incorporation  continues  in  the  bulk  of  the  film.  At  a 
nitridation  temperature  of  800‘  C,  the  peak  interfacial  nitrogen  concentration 
is  observed  to  remain  at  the  interface  for  nitridation  times  up  to  4  hours.  At  a 
nitridation  temperature  of  1150'  C,  the  peak  in  the  interfacial  nitrogen 
concentration  occurs  away  from  the  interface  for  a  nitridation  time  as  short  as 
5  minutes.  The  peak  interfacial  nitrogen  concentration  occurs  about  the  same 
distance  from  the  interface  as  is  observed  at  a  nitridation  temperature  of 
1000'  C,  but  saturates  at  a  lower  value. 

Modeling  of  the  kinetics  of  the  nitridation  has  been  done,  and  reveals  a 
mechanism  in  which  diffusing  species,  initially  consisting  primarily  of  the 
nitrogen  species,  react  with  the  substrate,  forming  a  nitrogen-rich  oxynitride  at 
the  interface.  At  sufficiently  high  nitridation  temperatures,  an  oxygen-rich 
oxynitride  subsequently  forms  at  the  interface  due  to  reaction  of  the  substrate 
with  an  increasing  concentration  of  diffusing  oxygen  which  has  been 
displaced  by  the  slower  nitridation  of  the  oxide.  The  data  are  consistent  with 
the  mechanism  discussed  above,  provided  the  influence  of  the  interfacial 
strain  on  the  nitridation  and  oxidation  kinetics  is  explicitly  accounted  for. 
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Specifically,  the  effects  of  interfacial  strain  must  be  taken  into  account  to 
explain  the  observed  width  of  the  peak  in  the  interfacial  nitrogen 
concentration.  This  is  an  example  of  a  chemical  inhomogeneity  being 
affected  by  a  structural  inhomogeneity.  The  fact  that  the  XPS  data 
qualitatively  confirm  the  predictions  of  the  kinetic  model  is  an  indication  that 
the  essential  physics  of  the  nitridation  process  is  contained  in  the  model, 
although  quantitative  discrepancies  indicate  that  the  kinetic  parameters  are 
not  yet  optimized.  In  spite  of  the  discrepancies,  the  model  successfully 
predicts  that  (a)  the  maximum  nitrogen  concentration  in  the  interfacial  region 
initially  occurs  at  the  interface  and  occurs  away  from  the  interface  with 
increasing  nitridation  time  at  sufficiently  high  nitridation  temperatures,  (b)  the 
maximum  nitrogen  concentration  in  the  interfacial  region  remains  at  the 
interface  even  for  long  nitridation  times  at  low  nitridation  temperatures,  and 
(c)  the  peak  interfacial  nitrogen  concentration  is  lower  for  the  highest 
nitridation  temperature  studied.  The  simulations  provide  valuable  insight  into 
the  reasons  for  these  experimentally  confirmed  observations,  and  also 
provide  insight  into  the  role  of  interfacial  strain  in  the  nitridation  kinetics. 

The  behavior  of  |AVFBj,  the  magnitude  of  the  flat  band  voltage  shift,  is 
correlated  to  the  observed  nitrogen  distributions.  This  correlation  is  shown  to 
consistently  follow  from  the  simple  postulate  that  it  is  the  strain-dependence 
of  the  kinetics  of  formation  of  defect  states  which  allow  lowering  of  the  system 
energy  by  relieving  interfacial  strain.  The  distribution  of  defects  is  deduced 
from  the  intensity  of  a  fluorine  marker,  and  is  found  to  be  consistent  with  the 
proposed  model  explaining  the  behavior  of  flatband  voltage  shifts  as  a 
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function  of  nitndation  conditions.  Macroscopic  intrinsic  stress  is  not  cbserveo 
to  affect  the  nitrogen  distribution,  but  does  affect  the  defect  formation  kinetics. 
Thus,  it  can  be  seen  that  the  observed  nitrogen  distributions,  the  kinetics  of 
the  nitndation,  and  the  electrical  properties  of  the  nitrided  oxides  can  all  be 
explained  in  a  self-consistent  manner  by  taking  into  account  known 
properties  of  the  initial  SiC^/Si  structure. 

Studies  of  SOI  structures  formed  by  ion  implantation  show  that  the  buriec 
cxioe, Si  interfaces  are  smoother  with  higher  implantation  temperature  in  the 
range  400*  C  to  525*  C.  While  the  stoichiometry  of  the  native  oxides  suggests 
that  a  higher  implantation  temperature  results  in  a  rougher  surface,  the  XPS 
data  on  chemical  depth  profiles  of  the  thermal  S^/surface  Si  interface 
snows  no  effect,  implying  no  effect  on  device  characteristics.  SOI  structures 
formed  by  zone-melting  recrystallization  have  also  been  studied,  but  these 
studies  have  not  yielded  useful  results  io  date. 
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Appendix  A 

Intensity  Analysis  of  XPS  Data 


In  the  present  work,  the  distribution  of  nitrogen  in  mtnded  oxide  films  s  me 
desired  auantin.  An  inhomogeneous  nitrogen  distribution  as  a  function  c< 
deDtn  m  the  film  is  evident  in  the  behavior  of  the  O  1  s  and  N  1  s  intensities," 
as  discussed  in  Chapter  3.  However,  this  behavior  is  biased  bv  me 
averaging  effect  of  the  relatively  long  photoelectron  mean  free  paths,  and  by 
the  fact  that  N  is  and  O  Is  photoelectrons  have  different  photoionization 
cross  sections,  and  have  signif:cantly  different  kinetic  energies,  and  hence 
different  mean  free  paths.  An  unbiased  measure  of  the  nitrogen  distribution  is 
required,  in  the  following,  the  simpler  case  of  a  homogeneous  film  will  first  be 
considered,  then  the  generalization  to  the  case  of  an  inhomogeneous  film 
composition  w ill  be  presented. 

The  intensity  of  photoelectrons  of  a  given  kinetic  energy  emitted  from  a 
homogenous  material  is  a  well  known  function. 17'25'34-36'1 17'124  Let  Ex  be  the 
x-ray  energy,  aN1sthe  N  Is  photoionization  cross  section,  DN  ox  the  atomic 
number  density  of  N  in  the  oxynitride  (SiOxNy),  and  ^nisox(En1s)  the  inelastic 
mean  free  path  of  N  1  s  photoelectrons  with  kinetic  energy  ENl5  in  SiOxNy.  In 
the  present  work,  Ex,  and  hence  the  kinetic  energy  of  photoelectrons  from  a 
given  core  level,  does  not  vary.  The  explicit  dependence  of  o  and  /.  on 
energy  will  therefore  no  longer  be  written,  but  rather  will  be  implicit  in  the  use 
of  different  k's  for  different  core  levels.  The  intensity  of  N  Is  photoelectrons 
from  SiOxNy,  dlN1sox,  which  originate  between  z  and  z  +  dz  below  the  surface 
is  given  by 
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dlN1s.ox  =  K^1SDN.oxexP[-z/^1s,oxSine]dz  ■  <A  '  ■ 

wnere  K  is  a  spectrometer  constant  depending  on  such  factors  as  x-ray  flux 
and  detector  efficiency  and  d  is  the  angle  between  the  plane  of  the  surface 
end  the  ^etectcr.  'rhe  dmect^d  intensity  of  ohotoelectrons  from  a 
semi-infinite  homogeneous  sample,  such  as  Si  2 p  photoelectrons  from  a  Si 
substrate,  would  be 


oo 

0 

* Si2p.Si  =  *Q  ^aSi2p^Si.SieXPf'Zy/^Si2p.SiSin®^dZ  =  ^°S.2p^Si.s/-Si2p.SiSin®-  *A2) 


The  total  detected  intensity  of  photoelectrons  from  a  film  of  finite  thickness 
such  as  N  1  s  photoelectrons  from  a  SiOxNy  film  of  thickness  tox  on  Si,  is 


^N1  s.ox^ox^  JQ  KctN1  sDN,oxexP['z/^Nt  s,oxS'ne]dz 
=  ^NtsDN.ox^NIs.oxS'^^-e^HoANU.ox3106])' 


(A3) 


The  total  intensity  of  the  substrate  Si  2 p  photoelectrons  detected  through  the 
SiOxNy  film  is 


L 15 


oo 


■sisas^ox)  =  J  KaS,2pDS<.S,exPHZ  •  t0x)^Si2p.S,S,ne]exP[-tox^S12p.oxS'ne]Cl2 


(A4) 


=  KaS,2DDSl.S1XS12aS,SineexP[-W>-Si20.oxSine]- 


Eq.  (A4)  can  be  divided  by  Eq.  (A2)  to  obtain 


,Si2p.Si(tox)/lSi2x».Si(°)  =  exPKA2p.oxSinel 


(A5 


For  notational  convenience  the  following  symbol  is  defined: 

AN  =  DN,ox^Nls.ox/DSi,Si^Si2p,Si  •  (A6 

« 

We  shall  hereafter  assume  that  the  intensities  are  all  divided  by  the 
appropriate  photoionization  cross  sections,  which  are  well  Known.31  Eq.  (A3) 
can  then  be  divided  by  Eq.  (A4)  to  obtain 

lNis.ox(tox)/ls,2p.s1(tox)  =  An  (1  -  exp[-tox/XNUoxsine])/exp(-tox/XSl2poxsine].  (A7) 

A  similar  equation  can  be  written  for  the  intensity  of  01  s  photoelectrons.  Note 
that,  from  Eq.  (A5),  tox  can  be  determined  if  ^N1sox  is  known,  or  XNlsox  can  be 
determined  if  tox  is  independently  measured.  DSi  Sl  and  ASl2pSi  have  13660 
measured, 46'49-75  so  that,  if  A.N1sox  is  known  or  can  be  estimated,  DNox  can  be 
determined  from  Eq.  (A7).  In  the  present  situation,  the  assumption  of  a 
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Homogeneous  film  is  invalid.  For  innomcgeneous  films  both  D  arc  /.  are.  m. 
general,  functions  of  z.  The  above  equations  must  therefore  be  generalized 
to  include  inhomogeneous  films.  We  begin  by  assuming  that  a  film  of 
thickness  tox  can  be  approximated  by  a  finite  number  of  layers,  each  of  which 
is  homogeneous.  The  i:h  layer  below  the  surface  can  then  be  characterized 
by  D  and  a,  (and  hence  A;).  and  by  its  boundaries  zh1  and  z,,  where  z0  (  =  0)  is 
the  surface,  cq.  (A3)  can  then  ue  rewritten  as 

>0, 

'nis.oA*)  =  K^DN.ox«Z>exP[-^N's.oAZ)Sine]C2 

"  f2' 

=  K  Z  D.nJ  exp[-zAiiN1SiOxsin0]dz 

1=1  z,-i 

n  (A8) 

=  Ksine  Z  Di,N,ox^i.N1  s.oxi^  '  exP[-(z,  -  Zi-l^i.NIs.oxS'H0]} 

i=1 

x  exp[-z,.rai  NlS0xs.n6]. 

The  thickness  of  the  i!h  layer  is  defined  by 


t,  =  z,  -z, 


i-i  • 


(A9) 


so  that 


n 

^OX  =  Xi  ^1  • 
1=1 


(A1 0) 
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The  z.„,  in  the  last  exponent  of  Eq.  (AS)  can  De  replaced  by 


i-1 


i-1 


2 


i-l 


(A1 1 


Equations  (A9)  and  (A1 1)  can  oe  substituted  into  Eq.  (A8)  and  t0  can  be 
detmed  to  be  0  to  give 


N1  s.ox'lox 


(to,)  =  Ksine 


I 

i=1 


^i.N.ox'S.NI  s.ox 


(i  -  expt-t/X; 


i.NI  s,ox 


sin©]) 


i-1  (A12) 

x  II  exp[-t:/Xj  N1  sin6] . 

J=0 


Similarly,  Eq.  (A4)  can  be  rewritten  as 


*Si2p,Si(tox)  “  KsinSDsj  si^-Si2p,Si  -TI  ©xP[“tj/Xj  Si2p.oxSin®l 


(A13) 


Equation  (A12)  can  be  divided  by  Eq.  (A13)  to  obtain 
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n 

^Nls,ox^ox^Si2p.Si^ox^  ~  N  ^  ’  ®*PHi'A-i  Nlsox3'^®]) 

1=1 

M  n 

X  n  exp[-t,X  N1S0xsine]  TI  exp[tk/Xk  Si2  xsin0]  (A14) 

j=0  '  •  k=0  ' 

Equation  (A14)  is  the  same  as  an  expression  previously  given36  (but  not 
derived),  but  generalized  to  include  the  energy  dependence  of  the  mean  free 
paths,  and  z  =  C  is  now  defined  at  the  film  surface,  rather  than  at  the 
film/substrate  interface  as  was  done  previously.36 

It  has  previously  been  shown119-121  that  the  mean  free  path  depends  on 
the  photoelectron  kinetic  energy  via  the  relation  X  a  ED,  where  p  is  in  the 
range  0.70  -  0.81  for  Si  and  Si02-  For  this  work,  p  =  0.75  has  been  chosen. 
In  Eq.  (A  14),  \N1s  ox  can  then  be  replaced  by  utilizing  the  relation 

\N1s.ox  =  ^N1?/^Si2p)°'75^i,Si2p,ox  ■  (A15 

For  notational  convenience,  the  following  symbol  is  defined: 

AN  =  (ESi2pENis)°'75-  (A1 6- 

Substituting  Eqs.  (A  1 5)  and  (A16)  into  Eq.  (A14),  the  following  relation  is 
obtained 
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^N1s.ox^oxK^Si2p.Si^ox)  “  2~*  -\n^  ‘  exPHiA|Ai,Si2p,oxSin®]) 

(> 

i-i  n 

x  I!  exp[-t,(AN  -  1)/Aj  Sl2poxsin6]  11  exp[V\s,2poxsine]  . 
j=u  k=i 

Equation  (At  7)  can  be  rewritten  as 

s.ox ^ox^^Si2p.Si^ox)  ~  ^Nts.ox^ox  *  ^Si2p,Si^ox  *  M 

(/ 

n 

+  . N ( 1  '  exp[-ti AN/X1  Si20.oxS|ne])  FI  exP[tA,Si2poxsinel  • 

i=1 

Equation  (A  1 3)  can  also  be  written  as 

lSi2p.Si^ox)/,|Si2p,Si^tox  -  M  =  exPHA, .Si2p.oxSin0]  ■  ( 

and  also  as 

n 

’s;2p.Si^ox)^Si2p.Si^^  =  n  exP['tAi.Si2p.oxSln®l  •  (■ 

Substituting  Eqs.  (A19)  and  (A20)  into  Eq.  (A18),  the  following  expression 
is  obtained 


\  _  ! )  /♦  \i  /♦  t  .  _  |  i*  ♦  \ |  /+  \\ 

~1'  ,N  _  •  N :  s  ox '  ‘ox ' 1  S>2p.S>  ’•  ’ov  ”  Nis.oxllcx  ■  Si2p.Si  lox'- 

an 

X  0si2f>.Si^),Si2f.Si(tox  "  '  (lS12p,Sl(,ox)/lSi2p.s/tox  *  M)  1} 


;A2‘ 


Ecuation  (A21 )  reduces  to  the  expression  previously  derived36  if  AN  =  1 .  i.e  r 
the  Kinetic  energy  deoenaence  of  the  photoelectron  inelastic  mean  free  oath 
can  be  negiectea.  Using  the  O  Is  photoelectron  intensities  anc  kinetic 
energy,  expressions  similar  to  Eqs.  (A15)  and  (A21)  can  be  obtained  for  A- 
and  A.  q.  It  is  to  be  noted  that  ihe  parameters  A,  N  and  A,  0  can  be  calculated 
entirely  from  the  known  quantities  AN,  and  A0  and  from  the  photoeiectrcn 
intensities  measured  at  film  thicknesses  of  zero,  tox,  and  tox  -  tv  A.  N  and  A.  D 
are  measures  of  the  composition  of  the  surface  layer  of  thickness  t..  Tne 
relative  concentration  of  N  in  the  surface  layer  is  directly  given  by 

WHO]  *  [ND  -  D.  N  o/(D,  Ooi  *  D,  N  01)  .  ANA.  n/(A0A,  0  *  AnA,  n)  (A22 

Measurement  of  the  photoeiectron  intensities  for  n  different  values  of  zr 
including  Zj  =  tox  (i.e.  the  substrate  intensity  with  no  overlayer),  yields  n  -  1 
values  of  A,  N  and  A,  0,  and  thus  the  relative  concentration  of  N  as  a  function  of 
depth  in  the  mtrided  oxide  film  can  be  obtained.  It  is  to  be  noted  that  a 
measure  of  the  homogeneity  of  the  film  can  also  be  obtained  at  a  fixed 
thickness  by  varying  0.30'34  or  by  varying  the  x-ray  photon  energy.29  To 
obtain  a  thickness  scale,  note  that  Eq.  (A20)  can  be  rewritten  as 
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V\si2p.ox  -  *Sin®  ^n^Si2p.Si^  fi^S^p.Si^^ 
1=1  1=1 


To  obtain  the  depth  scales  shown  in  Figs.  3.1  tc  3.8.  it  was  assumed  that 
XSi2oox  was  not  a  function  of  depth  in  the  film,  and  that  XSi2po)(  =  ^s.2p3C2- 
The  values  reported  for  X s,2psi02  w(ien  =  f  4-86.6  eV  are  generally  in  the 
range29’36'"22  27  to  30  A,  although  values  as  high  as  36  A  have  been 
reported29’"23  for  some  S1O2  films.  For  this  work,  ^-s,2pox  =  29  A  was  used  tor 
the  purpose  of  obtaining  the  depth  scales  in  Figs.  3.1  to  3.8  via  Eq.  (A23). 
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